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^\ • Abstract 

We extend the results and methods of Q to prove the existence of constant positive scalar 
curvature metrics g which are complete and conformal to the standard metric on S N \ A, 
where A is a disjoint union of submanifolds of dimensions between and (N — 2)/2. The 
existence of solutions with isolated singularities occupies the majority of the paper; their 
existence was previously established by Schoen jl^] , but the proof we give here, based on the 
techniques of is more direct, and provides more information about their geometry. When 
A is discrete we also establish that these solutions are smooth points in the moduli spaces of 
all such solutions introduced and studied in jjj and Q 

T^j- ! 1 Introduction and statement of the results 
O 

In this paper we construct solutions of the Yamabe problem on the sphere (S , go) with its 
standard metric which are singular at a specified closed set A c S N . More specifically, we seek 
a new metric g conformal to go and complete on A C S N , and with constant scalar curvature 
R. The problem may be translated into a differential equation as follows. Since g is conformal 

4 

to go, we may write g = u N - 2 go where u is a positive function on M \ A. The scalar curvature 
^Jy functions R{go) of go and R(g) of g are related by the equation 



On 



X 



4(N - 1) " u/ 4(iV - 1) 

In order that g be complete on S N \ A it is necessary for u to tend to infinity rather strongly on 
approach to A, and of course, we wish to solve this equation with R a (prescribed) constant. 
■ The first two terms of the operator on the left in (]l|) , namely 

N - 2 

Am = A SD - A{N _ l) R (9o) (2) 

give a second order linear elliptic differential operator known as the conformal Laplacian of the 
metric go- It satisfies the conformal equivariance property that if two metrics are conformally 
related, such as g and go above, then for any function </>, 

C go {u<t>)=u^L g {4>). (3) 

Notice that ([!]) corresponds to the special case of (§) when 0=1. 

This 'singular Yamabe problem' has been extensively studied in recent years, also in the case 
when the ambient manifold is more general than the sphere, and many existence results as well 
as obstructions to existence are known. Briefly, for a solution to exist on a general compact 
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Riemannian manifold [M, go), the size of A and the sign of R must be related to one another : 
if a solution exists with R < 0, then dim(A) > while if a solution exists with R > 0, then 

dim(A) < N ~ 2 and in addition the first eigenvalue of the conformal Laplacian of go must be 
nonncgative. Here the dimension is to be interpreted as Hausdorff dimension. Unfortunately, 
only partial converses to these statements are known. In particular, almost all of the existence 
results require A to be a submanifold, or at worst a stratified set. 

The character of the analysis required to prove existence of solutions when R < is quite 
different than in the other cases. In fact, it is not too difficult, using barrier methods, to construct 
solutions of The difficulty is in showing that these solutions have a singularity at A strong 
enough for g to be complete on the complement. Similarly, when R = 0, ([j]) becomes linear, 
albeit of a somewhat nonstandard sort. We shall not discuss these cases any further, but refer to 
the work of R. McOwen |[o]], D. Finn [Q, |2j and the first author ||, and references therein. 

We shall concentrate here on the case R > 0. The first examples of solutions were constructed 
by R. Schoen |l2| when A is either a finite set of points (of cardinality greater than one) or a 
nonrectifiable set arising essentially as the limit set of a Kleinian group. Nontrivial solutions 
with higher dimensional singular sets were constructed by the first author and N. Smale Q, the 
second author [jll), and finally in some generality by both of us ||. This last result states that 
solutions may be constructed on an arbitrary compact manifold of nonnegative scalar curvature 
M whenever A is a finite disjoint union of submanifolds of dimensions between 1 and T 2 . This 
paper is meant to extend the methods of Q to also allow for the construction of solutions with 
isolated singularities in the case where M = S N and go is the standard metric on S , i.e. to 
allow A to be an arbitrary disjoint finite union of submanifolds of dimensions between and N ^ 2 . 
Our main result is: 

Theorem 1 Suppose that A is a disjoint union of submanifolds in S N of varying dimensions. 
Write A = A' U A", where A' is the union of all submanifolds of dimension zero, i.e. A' = 
{pi, . . . ,p n } is a collection of points, and Lambda" = U" =1 A? w ^ ere dimAj = kj. Suppose 
further that < kj < N ^ 2 for each j , and either n = or n > 2. Then there exists a complete 
metric g — u N - 2 go on S N \ A which has constant positive scalar curvature R = N(N — 1). 

For most of this paper we shall consider the case where A" is empty, so that A is a finite 
collection of points 

A = { Pl ,..., Pn }cS N . 

The modifications needed to treat the general case are an amalgamation of the techniques here 
and those of [||, and shall be described in a brief final section of this paper. 

Of course, solutions with isolated singularities were already constructed by R. Schoen, but 
his remarkable proof is long and difficult. We feel that (the extension of) the methods of 
avoid many of the difficulties he encountered and substantially clarify the analysis. In addition, 
various properties of the solutions follow immediately from the construction here, but are not at 
all obvious for Schoen's solutions; we describe some of these below. 

Let us describe the two approaches, and some of their differences, somewhat more specifically. 
If we normalize the desired constant scalar curvature R to be N(N — 1), then the equation we 
are trying to solve on the sphere is 

.r, s . N(N-2) N(N-2) n+2 

Af(u)=Au y — J -u+— y — '-u— =0. (4) 

In order to construct a solution with a given singular set A, first a one-parameter family of 
approximate solutions u e of (^), each element of which blows up quickly enough near A, is 
constructed. Then (^) is linearized about u £ , and the resulting linear operator L £ is analyzed. 
If it can be proved that L e is surjective on some reasonable space of functions, at least when 
e is sufficiently small, then a standard iteration argument may be used to correct u e by adding 
to it a function v to obtain an exact solution to (nT) which blows up sufficiently quickly at A. 
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Unfortunately, L e is not surjective on L 2 ; in fact, L £ is self-adjoint on L 2 , but is in its continuous 
spectrum, so that it does not even have closed range on this space. Schoen's tactic is to find an 
explicit infinite set of functions which span an approximate nullspace K, such that L e is invertible 
on K- 1 . He first solves the equation on K ± , and then gives 'balancing' conditions to ensure that 
the solution of this restricted problem is a solution of the original problem. Our somewhat 
different tactic is to work on a finite dimensional extension of a certain weighted Holder space. 
On this space, L e is actually an isomorphism, when A is discrete, and surjective in general, for e 
sufficiently small. Unfortunately, as e tends to zero, the norm of any right inverse for L e blows 
up (Schoen encounters a similar problem). We must then analyze precisely the rate and manner 
of blowup. This is different than in fj], when A' is empty, and where the (right) inverse for L £ 
is bounded as e — > 0. The advantage of working in weighted spaces over using the approximate 
nullspace K is that in our approach there are only a finite number of obstructions to solving the 
equation, and these may be identified explicitly and geometrically. In fact, these obstructions are 
intimately connected with the definition of the finite dimensional extension mentioned above. 

As noted earlier, another advantage of our approach is that we may easily derive various 
properties of the solutions. The main property we are interested in is nondegeneracy, which 
will be defined precisely in §11. This property is important in the study of the 'marked' and 
'unmarked' moduli spaces Ma and M n of solutions of this problem where the singular set A 
is fixed or allowed to vary amongst all configurations of n points in S N . There are notions of 
marked and unmarked nondegeneracy associated to each of these spaces. These moduli spaces 
were defined and studied by the first author, with D. Pollack and K. Uhlenbeck in {jj, |§. It 
was proved that for M = S N they are real analytic sets. If there exists some g € Ma which is 
marked nondegenerate, then the top dimensional stratum in the component of g is a real analytic 
manifold of dimension n. Similarly, if g £ M n is unmarked nondegenerate, then nearby g this 
moduli space is a real analytic manifold of dimension n(N + l). Nondegenerate solutions on S of 
a very special type were constructed in S ; for these solutions, only certain special configurations 
A (in particular, only those with cardinality n an even number and with points clustered in pairs) 
are allowed. In contrast, we prove 

Theorem 2 For any integer n > 2 and any configuration A of n points in S there exists an 
element g 6 M n which has singular set A and which is unmarked nondegenerate. For a generic 
(in fact, Zariski open) set of A, this solution is marked nondegenerate, and for such a metric the 
points (pi, . . . ,p n ) of A and the 'Delaunay necksizes' (ej., . . . , £„) constitute a full set of coordinates 
in M n near g, while the Delaunay parameters alone yield coordinates in Ma near g. 

The Delaunay necksizes will be defined in §2. We remark that the admissible sets of Delaunay 
necksize parameters {ei, . . . , s n } in the construction are not arbitrary sets of small numbers, but 
are required to satisfy a 'balancing conditions' (|l2]). 

Although the statement of Theorem [j] is for metrics on the complement in the sphere of a set 
A, we shall use the conformal equivariance of the equation (Q) and prove instead the existence of 
solutions of this equation on the complement in M. N of a (finite) set S which decay at infinity. 
This makes the technical details somewhat easier in that we always have a preferred coordinate 
system. 

The solutions we construct here are not 'the same' as the ones constructed by R. Schoen; they 
are also quite different from the ones constructed in |8| , and possibly don't even lie in the same 
components of the relevant moduli spaces. The ones in the latter paper are not required to have 
small Delaunay parameters, while (at least some of) the ones in jl2| may be thought of as infinite 
strings of spheres, connected by very small necks, and joined together at a central convex sphere. 
Our solutions also have small necks, but the central region is metrically both concave and very 
small. 

Karen Uhlenbeck has informed us that she too has established the existence of solutions to 
this problem with exactly one singular point on compact Riemannian manifolds (M,go) with 
nonncgative conformal Laplacian. The construction of the approximate solutions in this case is 
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substantially different than ours, and requires the positive mass theorem, but her linear analysis 
is in roughly the same spirit as ours. 

The connection between this problem and the construction of complete, noncompact surfaces 
in R 3 with constant mean curvature (CMC) is well-known. Indeed, Kapouleas' initial construction 
of these surfaces || is related to Schoen's construction of constant scalar curvature metrics fl2|| , 
and the analysis of the moduli space problem in jl| for the CMC case was directly inspired 
by 0. Maintaining this tradition, we shall show, in a work in preparation, the existence of 
complete, noncompact CMC surfaces in R 3 with an arbitrary number of ends (greater than two), 
and with prescribed Delaunay parameters on these ends. These solutions are nondegenerate, 
hence are smooth points in the relevant moduli spaces. The geometric difference between these 
solutions and those of Kapouleas are the concavity in the central region of our solutions, versus 
the convexity of his. 

The outline of this paper is as follows. In the first section we examine the basic models of 
singular solutions, the Delaunay solutions (the reason for this name is explained in Q), and prove 
some estimates on them that will be required later. After this, we present the construction of 
approximate solutions for the problem. These functions are periodic near each end in a certain 
natural coordinate system, and correspondingly the linearized operators we must analyze have 
periodic coefficients. Although an analysis of such operators was made in JtJ, we need rather 
stronger behaviour than seems to be immediately available from the methods there, so the next 
few sections are devoted to a somewhat novel construction of inverses for the linearized operators. 
First of all, the linearized operators are invertible only on rather special function spaces. Then, we 
construct an inverse for the Dirichlet problems for the linearizations in the e-neighbourhoods of 
the singular points, and on the complement of these neighbourhoods. By analyzing the Dirichlet- 
to-Neumann operators on the union of the boundaries of these neighbourhoods, we may join these 
right inverses together to get a global right inverse. To carry this out, we take a slightly indirect 
path and prove the existence of inverses by this method only for a model problem; the existence 
of inverses for the true problem is deduced by a rather intricate sequence of perturbations. Af- 
ter proving sufficiently fine estimates for these inverses, a rather standard contraction mapping 
argument (estimates for which, unfortunately, are not so trivial) is employed to complete the 
construction. After this we discuss the issue of nondegeneracy of the solutions. As noted earlier, 
in most of the paper we discuss only the case where A is discrete, but in the last section we discuss 
the changes needed to handle the general case. 



2 Delaunay type solutions 

In this section we recall some well known fact about the Delaunay solutions of ([!]) on S N that 
will be used extensively in the rest of the paper. A reference for facts not proved here is 0, 
cf. also It is known that if w is any solution of (Q) on S N \ {pi,p2}, then it is invariant 

with respect to any conformal transformation fixing these points; if these points are antipodal, 
as may be assumed, then u is rotationally invariant. In either case, the equation it satisfies may 
be reduced to an ODE. It is convenient to stereographically project the sphere S N to M. N from 
one of the singular points, say pi, so that P2 is sent to € Mr. Then the solution is transformed 
to a radial solution of 

,/V N A N(N-2) N + 2 

Af(u) = Au + — — '-u— = 0, (5) 

on 1^ \ {0} which is singular at the origin. 
We may reduce (^) by writing 

2 — N 

u(x) = \x\^~v(—\og\x\), 
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and using t = — log x to get 

d 2 v (N-2) 2 N(N~2) n± 



-1! - 



V 7T=2 = 0. (6) 



dt 2 4 4 

This equation is nondissipative, and the Hamiltonian energy 

„, . 2 (iV-2) 2 2 (iV-2) 2 j2N_ 

H(v,v)=v 2 -- —v 2 + - J — v «-2 (7) 

4 4 

is constant along solutions of (Q). Here and in all that follows, • denotes differentiation with 
respect to t. By examining the level curves of H, we see that all bounded solutions of (0) lie in 
the region of the (v, u)-plane where H(v,v) < 0. There are several different types of bounded 
solutions; we summarize their basic properties: 

Proposition 1 For any Hq G (— ^y^{^j^) N / 2 , 0), there exists a unique bounded solution of (Qj 
satisfying H(v,v) = Hq, v(0) = and v(0) > 0. T/iis solution is periodic and for all t £ 1 we 
/iave u(t) G (0, 1). TTiis solution can be indexed by the parameter e — v(0) G (0, (-^r^)^ -2 -^ 4 ) , 
which is the smaller of the two values v assumes when i) = 0. When Hq = — jV 2 ~ 2 ( ) N / 2 , £/iere 
is a unique bounded solution of given by 

2 — JV 

Finally, if v is a solution with Hq = i/ierc either v(t) = or u(t) = (cosh(i — ^o)) - ^ /» r some 

< e E. 

We will call these the Delaunay type solutions. Although we do not know them explicitly, the 
next two propositions give sufficient information about their behaviour as e tends to zero for our 
purposes. 

Proposition 2 Fix e G (0, ( JV ~ 2 ^)( jy ~ 2 )/ 4 ) and let v £ be the corresponding Delaunay solution. 
Then the period T £ of v e tends to infinity monotonically as e — > and satisfies 

Te = - (j^ + log(e). 

In addition, for all t El 

, ,N-2 s 
e < u e(^J S ecosh( — - — t) 



Proof : The second claim is rather simple. Since H is constant along solutions, 

„„,„„,, . a - m„, ; _„«, . _fi^2£y - ,8) 

Moreover, it is easy to see that w e is increasing for t G [0, T e /2] and decreasing for t G [T e /2,TJ. 
In particular the minimum of is achieved at t = and equals e. Using this lower bound in (g) 
we see that 

. 2 (iV- 2) 2 2 2^ / ^^/(^-2) 2 , 2 2n 

Taking the (positive) square root and integrating this differential inequality yields the inequality 
v £ (t) < ecosh(^t). 
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Now we derive the asymptotic estimate for the period T s . Let w ma x,e denote the larger of the 
two solutions of 



9 ±±L 9 ^iv 

v - v N - 2 = e — e N - 2 . 



(Hence w max ,e is the larger of the two values v e assumes when v s = 0) . By M) , we get 



T e /2 



2N 

N=2 _2 



s ~ + eT^y 1/2 v E dt 



N -2, 



In the interval of integration we may change the variable of integration to v = v £ (t)/e. This yields 



tWx.s/e 4 j •. wo, -/V — 2 m 
(v 2 - 1 - _ l))- 1 /^ = T e 



It is then not difficult to see that this integral grows as — (1 + o(l)) loge. A proof of the fact that 
the period is a monotone function of e is indicated in M and also proved by elementary ODE 
methods in an unpublished manuscript of Chouikha and Weissler. □ 

Proposition 3 For any e 6 (0, (^j^ 2 -)^^ 2 ^ 4 ) and /or any f G K the Delaunay solution v e 
satisfies the estimates 



v e (t) — e cosh 



N-2 



<cjv£ N - 2 e 2 ' 



. . TV — 2 , ( N — 2 \ «+2 n+2, 

u e (t) — esinhl — - — il <c N e Tr ^e 2 1 

/or some constant C/v > which depends only on N . 
Proof : We start with the a priori estimate 

W-W-2 , JV-2 , 

u e (*) < ecosh( — - — t) < ee 2 , 

valid for t > 0. Next, writing the equation for v £ as 

d 2 v (N-2) 2 N(N-2) n+2 



df 2 



we can represent v e as 



N - 2 iV(7V - 2) jv 



u £ (i) = ecosh( — - — t) 

This leads immediately to the estimate 

,N -2 



e 2 



,(2-N)s 



N-2 JV + 2 

e^~ z v £ N ~ 2 (z) dz ds. 



< ecosh( — - — t) - v e (t) < 



N — 2 N + 2 N + 2 



■ N-2 p 2 



Finally, differentiating the formula of v £ with respect to t, we get 



, , . N-2 . /N-2 
v e(t) 7. — esmh | — - — t 



< 



N(N-2) (N-2 N-2 t 



,(2-N)s 



N-2 M+l 

e— z v £ N - 2 (z)dzds 



+e~ 



N-2 N + 2 



e 2 z v n-2 (z)dzds < 



N 2 -A 



IG 



N + 2 N + 2_ + 
- £ N-2p 2 1 



□ 
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Remark 1 The estimate of this Lemma is only interesting in the domain where 

N + 2 N + 2 I .1 , N— 2 u| 

2 CAr£ iv^2 e — 1*1 < (JV- 2)ee— , 
i.e. where \t\ < — N 2 _ 2 loge — c^v, which is close to half a period of v s . 

We transfer these estimates and this remark back to the x-coordinates to obtain 

Corollary 1 For any e G (0, '' ] and any x G K \ {0}, the Delaunay solution u e {x) 

satisfies the estimates 

F . „ ... N + 2 , 

< C N e~^\x\ 



rd r u £ (x) + ^-l £ \x\ 2 - N 



< C N e N ~ 2 \x\ 



-N 



In particular, in the annulus e N - 2 < \x\ < 1, u e is well-approximated by 



12 — JV 



There are some important variations of these solutions, leading to a 2/V+2-dimensional family 
of Delaunay type solutions. These variations are families of solutions U(s) of Af(U(s)) = with 
U(0) = u e , depending smoothly on the parameter s. The derivatives of these families with respect 
to s at s — correspond to Jacobi fields, i.e. solutions of the linearization of JV about one of the 
u £ , and will be described below. 

We describe these families of variations in turn. The first is quite obvious. It is the family 
where the Delaunay parameter e is varied 

(-£, 1 - e) 3 ?7 — > u s+v (x). 

2 — N 

The second corresponds to the fact that if u is any solutions of J\f(u) = 0, then R~^~u{x/R) also 
solves this equation. Applying this to a Delaunay solution yields the family 

R+ 3 R — > |x|^w £ (-log|x| +logi?). 

The other two families of solutions correspond to translations. The simpler of these is the usual 
translation 

R N 3 b -> u £ (x + b). 

The final one corresponds to translations 'at infinity'. To describe this we use the Kelvin trans- 
form, u(x) — > \x\ 2 ~ N u(j^2 ), which preserves the property of being a weak solution of (||). Start 
with a Delaunay solution u e (x). Its Kelvin transform is 

\x\^-v e (log\x\). 

Translate this by some a G M. N to get 

\x - a | (log | a: - a\), 
which has its singularity at a rather than 0. Its Kelvin transform yields the family 

R N 3 a -> \x - a|x| 2 |^ u £ (-21og|xj +log|a; - a|a;| 2 |). 



Each function in this family has singularity at again. 

In fact, the first and third variations, where the Delaunay parameter is changed or where the 
singularity is moved, are less well behaved than the other two: the Jacobi fields corresponding to 
them grow too quickly. Thus we shall focus attention on the smaller family of solutions 



2-N 

1 2 I 



u e (R,a,x) = \x - a\x\ 2 \ 2 w e (-21og |x| + log \x - a\x\ A \ + logi?), (9) 
where only translations along the Delaunay axis and of the 'point at infinity' are allowed. 
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3 Construction of the approximate solutions 



We now proceed to define a family of approximate solutions to the problem, using the Delaunay 
solutions of the last section. Although we ultimately wish to construct solutions on S N \ A, we 
shall use the conformal invariance of the problem and construct solutions on ~R N \ E, where 

E = \X\i . . . , x n } 

with appropriate regularity at infinity. This is done purely to make the notation as simple as 
possible. 

The approximate solutions we define here have Delaunay type singularities in a ball B(xi,pi) 
around each point of E and are harmonic outside the union of slightly larger balls. We use not 
just the radial Delaunay solutions, but elements of the family u e (R,a,x) defined above. The 
additional parameters in this family yield, on the linear level, the Jacobi fields constituting a 
'deficiency space' needed to obtain surjectivity of the linearized operator. In addition, these extra 
parameters are essential in ensuring a sufficiently good agreement in the transition regions around 
the boundaries of the small balls. We emphasize again that we do not use all possible Delaunay 
solutions in each ball, but just two of the four possible families. 

Let 

e = (Ei,...,£„), R = (Ri, . . . , R n ), p=(pi,...,p n ), and a = (ai, . . . , a n ) 

be sets of (small) Delaunay parameters, positive real numbers corresponding to translations of 
the Delaunay solutions, (small) radii and vectors in Mr, respectively. We shall impose various 
constraints on these quantities. First, assume that the balls B(xi, pi) do not intersect one another. 
We will denote by Sip the set R N \ U2 =1 B(xi, pi). Next, we shall assume that 



Pi = e? • (10) 

Finally, fix a set of positive numbers q±, . . . , q n , and assume that 

Ei=eqi, i=l,...,n, (11) 

or simply s — eq, for some e > 0. We shall call a set e satisfying ( [Tl| ) admissible. Notice that ( fl0| ) 
and ( |TT| ) imply that the quantities £i, . . . , £„ are all comparable to one another, and similarly for 
the pi. Because of this, it makes sense to estimate various quantities later by e or p = e 4 ^ N ~ 4 \ 
We shall now define the family of approximate solutions Ug(R,a,-). In each ball B(xi,pi), 
Ue(R,a,x) is equal to u £i (Ri,cn,x — Xi). On the exterior region Q,p, Ue(R,a,x) will equal some 
harmonic function w(x). With this definition, it is possible to make ug(R, a, x) continuous across 
each dB{xi, pi), but not C 1 . For example, if one takes w(x) to be the unique harmonic function 
which decays at infinity and whose value on dB(x i: pi) is given by u £i (Ri i) then, by 

definition u s (R,a,x) is continuous but not C . It is easier to deal with a harmonic function w 
which is a sum of multiples of Green functions with poles at the Xi. The next result shows that 
we can mitigate this discontinuity of derivatives, at least to the extent needed later, by choosing 
all parameters carefully. 

Proposition 4 Suppose that e is a set of admissible Delaunay parameters, with e < Eq. Suppose 
that the parameters Ri and ai satisfy the relationships 

Y^RT^Rff^-x^^q^ i = l,...,n, (12) 

and 

1 N ~ 2 \ y W-2 

a *o = - R t 2 2-^liRz 2 \Xi a ~Xi\~ (x io -Xi), (13) 
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and that the harmonic function w s (R,x) on fip is given by 



w g {R, x)=Y. i R ^~ \ x ~ X ^ N - ( 14 ) 



2 

i=i 



Then, for i = 1, . . . , n and j = 0, . . . , N, we have the estimates 

(u Si (Ri,ai,- -Xi) - w s (R,-))(xi + pi9)4>j{9) d9 = 0(ep 2 ), 



and 

Here 4>j(9) — 9j, j — 0, . . . ,N, are eigenfunctions corresponding to the first two eigenvalues of 
the Laplacian on S l , i.e. Xq = and Ai = N — 1. 

Remark 2 The definition ( flj| j says that q is the eigenvector of some matrix, with coefficients 
depending on the Ri and Xi, with eigenvalue 1. Thus we first fix q, and then only use parameters 
Ri for which ( jll j is valid. Appropriate choices of Ri do not exist for every q, but it is not hard 
to show that there is a cone in the positive orthant of the q space for which solutions exist. The 
quantities R and a will be considered fixed and bounded (with each R^ 1 also bounded), so we 
usually neglect to mention their appearance in various constants in the estimates below. 

Proof: To prove these estimates, we shall expand both u ei (Ri, ai, x — Xj) and vj s (R,x) near 
dB(xi,pi). First we consider u ei (Ri i). From Proposition^, if \a\p < 1/8 as we may well 

assume by taking p small, then near dB(xi, pi) we have 

£ ■ / N-2 2~N \ 

u ei (Ri, ai ,x-Xi) = ^[R i 2 \ x - Xi \2-" + R. 2 (l + (JV-2)(ai • (x-x t ))J + 0(ep 2 ). 

It is precisely at this step, in estimating the size of the error term, that we are led to fix the size 
of the balls using ([To|), since with this definition of p we have ep 2 = e N - 2 p~ N ■ We also have 

N — 2 JV-2 2-iV 

d n u ei (Ri,ai,x - Xi) = — e i (R i 2 \x - Xi\ ~ -R t 2 \x - Xi\~ (a» • (x ~ Xi))) + 0(ep). 

Next, for the expansion of w s (R, x) near dB(xi , pi ) we have 

w e -(R, x) = £ -fR7f \x - x l0 \ 2 - N + l^fnizio " Xi\ 2 ~ N + 



and 



(2 - A0|x lo " x i\ N ( x io - x i) ■ ( x ~ x io)) + °( £ /° 2 )> 
d ri w s (R,x) = (2- N)(^-Rff\x - x t0 \ 1 - N + 
^ ^R^~\xi - x;| _JV |x - x <0 | _1 (x <0 - Xi) ■ (x - Xi,,)) + 0(ep). 
Taking the difference of these expressions we obtain, near |x — Xj | = pi , that 



£' 2 — — v — > £■ N-2 

Ue io {Ri ,a io ,x-Xi ) -w s {R,x) = I -yi? l0 2 - 2^ -j-Ri 2 \x l0 -Xi\ 2 ~ • 
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+ (AT - 2) [ £ -fR*fa i0 + ]T ^R^\xi 



Xi\ N {x io -Xi) I • {x - x io ) + 0{ep 1 ) (15) 



and 



9 n a (u £io (Ri ,a io ,x-x io ) -We(R,x)) =(N-2)[ ^-R H ? a lQ + 



(16) 



-AT, 



|x 1 



Now insert jl^ ) and (^) into these expressions, using ([n]) again, to obtain the result. □ 
Finally, using the function Wg(R, x) given in Proposition 0, we define the approximate solution 



i(R,a,x) — ^2 Xi( x ~ Xi)u £i (Ri,ai, x - xf) + w g (R,x)(l - ^ \i(x - Xj)). 



(17) 



Here the \i are smooth, radial functions equalling one in \x\ < pi, vanishing in \x\ > 2pi, and 
satisfy \d r Xi(x)\ < cp^ 1 and |c^Xj(x)| < C P7 2 ■ Notice that we may define the approximate 
solution Ug(R, a, x) regardless of whether the relations (|l^) and ( |l3| ) are satisfied. 

We conclude this section by stating an estimate which will be used extensively later, and 
which follows easily from the proof of Proposition ^. 

Corollary 2 Suppose either that a satisfies J7j[ ) or a — 0, and that R satisfies Then for 

Pi < \x — Xi\ < 1 



u g (R, a,x) N - 2 - u ei (Ri,a,i,x - Xi) N - 2 < Ce N - 2 (p\x - x { \ N 6 + \x~x l \ N b ) 



|JV — 5\ 



Proof: Write Ug and u £i for the functions appearing in the estimate. Also, assume x% = for 
convenience. We write the quantity to be estimated as 

4 4 _ _ 4 i_ 

\(XiU Ez + (1 -Xi)wg) N - 2 -uil 2 \ = \(wg + Xi(u Ez -wg))"-z -nil 2 |. 
By Taylor's theorem, this is dominated by 

C(tUg)^ _1 |u e » - w g \. 

From (|l5|), we see that when R solves (|l2|), then \u £i — Wg\ is bounded by Cep 2 when a solves 
(|l3|), and by Ce(p + \x\) when a = 0. Now use that Wg < Cs\x\ 2 ~ N to complete the proof. □ 



4 The linearized operator on R \ {0} 

Much of the analysis in this paper concerns the linearization of the nonlinear operator Af about one 
of the approximate solutions Ug(R, a, x). We shall approach the study of this linearization grad- 
ually. In this section we define the linearization about one of the Delaunay solutions u £ (R, a, x) 
on WL N \ {0}. The main point here is the introduction of the Jacobi fields ^ J E '% a (x). In succeed- 
ing sections we discuss the refined mapping properties of the linearization, first for the Dirichlct 
problem on the unit ball, then for a simpler, locally radial model for the true linearization on 
M. N \ E, and finally for the true linearization itself. 

Fix one of the Delaunay solutions u £ ^ a - The linearization of Af at u £ ^ a is defined by 



T d 

as 



Af(u £ (R, a, •) + sv) = Av+ N( " N + 2 \ £ {R, a, -)t^v. (18) 



s=0 
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More generally, this operator can also be defined as the derivative at s = of J\f(U(s)), where 
U(s) is any one-parameter family of solutions with U(0) = u e (R,a,x), U'(0) = v. Viewed this 
way, it is immediate that varying the parameters in any one of the families of Delaunay solutions 
leads to solutions of L £ ji a ^ — 0. Solutions of this homogeneous problem are called Jacobi fields. 

The Jacobi fields corresponding to these various families are easy to determine rather explic- 
itly, cf. @. We shall be primarily interested, at least initially, in the Jacobi fields at a = 0. Since 
the functions u £ (R, 0, x) are radial, we may separate variables in this case. Introduce the eigen- 
functions <fij(0) and eigenvalues Xj for the Laplacian on S ,JV_1 . Then corresponding to Ao = are 
the Jacobi fields for the first and second types of variations: 

. *°-(a:) = ^(x) = ^^^(-loglxj) = |x|^$°-(-log|x|), 

. = -\x\^(x) - ^u £ (x) = N ^^(_log|*|) ee \x\^*° e ' + (-tog\x\), 

Here we take (j>o{&) = 1. The third and fourth families have Jacobi fields corresponding to the 
next set of eigenvalues, Ai = . . . = Xn = N — 1. These are: 

. = lar^C^l^^C-logla;!) - log ^(») = N^^'+C-log |o;|) ^-(fl), 

. n>-(x) = = \x\~f (l^v E (-log\x\) - ^(-log|x|))^(<?) 

= M^*i.-(-iog|s|)&(0). 

Later we shall also use the Jacobi fields corresponding to differentiating the family u £ {R, a, x) 
and evaluating not necessarily at R = 1, a = 0. We shall denote these Jacobi fields by 

*i'X» = \*\™H'Xa( X )> j = 0,-" ,N. 



5 Function spaces 

We shall consider the action of the linearized operator on weighted Holder spaces, which we now 
define. Since we shall be considering both local and global versions of this linearization, we first 
define the appropriate spaces on the ball: 

Definition 1 For each k EN, < a < 1 and a 6 R + set 

, IV7 , , J k+a f\W k u(x)-V k u(y)\\ 

\lt\ k ,a,ltr,2*] = Sup \ \+0- L+a Sup M \-t ■ 

M6M»] lpj J \x\,\y\ela,2a]\ F - V\ a J 

Then, for any |t £ R, the space C k,a (B(0, R) \ {0}) is the collection of functions u which are 
locally in C k,a and for which the norm 

\\u\\k,a,p = Sup cr _Al |M|fc. Q ,[ CT ,2 CT ], 
<y<R/2 

is finite. The subspace C'%(B(0 t R) \ {0}) consists of those functions which, in addition, vanish 
at the boundary \x\ = R. 

In particular, the function r M is in C k,a for any k, a or /i. 

There are analogous spaces of functions on the complement of the prescribed singular set in 
S N of the eventual solution defined in the obvious way. Since it is easier to work in Euclidean 
coordinates, we shall usually consider instead the induced operator and problem on M. N \ E, 
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where S = {xi, . . . , n}. The only difference is that we now need to consider decay properties near 
Euclidean infinity. Consideration of these is included in the following definition. Divide R N \ E 
into the union of three slightly overlapping open regions: ft± = U 7 j l =1 B(xi, 1), a neighbourhood of 
E, = {\x\ > C}, a neighbourhood of infinity, and a bounded piece f^. 

Definition 2 For any v, v' £ R, the space C^'", (M. N \ E) is defined to be the collection of u £ 
C k - a (R N \ E) for which the norm 

\\u\\k,a, u y = ||u||jfe,a,i/,ni + IIM _Ar+1/ u (x/\x\ 2 )\\k, a ,B(o,i/c) + \\u\\k,a,n 3 

is finite. 

In this definition \f\k, a ,n denotes the usual Holder norm of the function / over the (bounded) 
open set Q. 

The one result about these we shall use frequently, and without comment, is that to check if 
a function u is an element of some C°' Q , say, it is sufficient to check that d(x)~ v \u{x)\ < C and 
d(x)~ u ~ 1 \'Vu(x)\ < C. Here d(x) is a smooth positive function on the complement of E, equaling 
\x — Xi\ on each B(xi, 1). 



6 The linearized operator on the unit ball 

In this section we shall study the Dirichlet problem in the unit ball B(0, 1) for the linearization 
about one of the radial Delaunay solutions u £ (R, 0, x), i.e. when a — 0. We restrict to the radial 
Delaunay solutions because they may be studied using separation of variables. Later we shall use 
a perturbation argument to treat the more general case. 

For simplicity, denote the linearization Le^o by L £i r, and when R = 1, simply by L £ . Our 
goal here is to study the problem 

L s>R w = f in B(0,1)\{0} 

(19) 

w = on dB(0,l). 

Specifically, we wish to find suitable spaces on which L £t R is surjective, and to find optimal 
estimates for the solution of (O) in these spaces, uniformly as e tends to zero. The first part of 
this was already considered inH. There it was proved that L £y R is surjective as a map between 
certain weighted Sobolev spaces, but the uniformity in e was not considered. For various technical 
reasons, it is more useful here to consider the action of L £ r on the family of weighted Holder 
spaces defined in the previous section. 

One can prove, analogous to the result of @, that 

L e , R : CjS(B(0, 1) \ {0}) — C°l Q 2 (i?(0, 1) \ {0}) (20) 

is Fredholm for all fx £ R, ^ ^ {f-fi^)}, where {^f(s)} is a discrete set, with each nf(s) 

determined by the ordinary differential equation induced by L £ ^ on the eigenspace of the 
Laplacian on S 1 ^^ 1 . The values of the inadmissible weights f-fi 5 ) ma y depend on s (although, 
somewhat miraculously, it turns out that those with j = 0, . . . , N do not). These values play the 
same role as the the indicial roots of a regular singular problem; indeed, they arise in the same 
manner as the rates of growth or decay of solutions of L £ ^w = 0. However, they are determined 
in a far less explicit manner, analytically rather than algebraically, and so are usually impossible 
to compute. Fortunately, we know the solutions of this equation, hence the rates of growth of 
these solutions, explicitly when j = 0, . . . , N, so we can determine, as in @], that 

for all M <^-, ^,eC^ for ah> < (21) 
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and for j = 1, . . . , N, 



^'X- eC M°' forall M <i^, for all M < (22) 

From this one sees that ^ = (2 - N)/2 and ^+ = (4- N)/2, fij = -N/2, j = 1,...,N. As we 
pointed out above, the values /i^, for j = 0, . . . , N, do not depend on e. 

It turns out that ( p0| ) is surjective when /x < (2 — N)/2, [i ^ [ij , and injective when /i > 
(2 — N)/2. The basic conflict is that, although this indicates that we should be working on spaces 
with weight less than (2 — N)/2, we are using L £yR to correct the approximate solution, which 
only grows like r^ 2 ~ N ^ 2 : and we obviously do not want the correction term to blow up faster 
than this. As proved in this may be remedied by working in a finite dimensional extension 
of C 2 ' a , where (2 — N)/2 < /i < (4 — N)/2 and the extension is the span of the Jacobi fields 

\D'°'^. The point is that although in general ( |l9| ) has no solution w 6 Cy Q when / £ C°lf 2 and 
/i > (2 — N)/2, there is a solution with the correct decay but with eigencomponent wq having 
the wrong Dirichlet data. To correct this, one may add on some multiple of one or the other of 
the Jacobi fields We shall prove something similar, but with two additional considerations. 

The first is that unlike the analogous situation in [|| , any choice of right inverse has norm blowing 
up as e tends to zero. However, this blow-up happens only on the subspace of functions spanned 
by the eigenfunctions <fij(0), j = 0, . . . , N. In order to obtain reasonable estimates, we need to 
work in the space with weight /i with 1 < fj, < 2. For this reason we shall also need to append to 
this 'deficiency space' the Jacobi fields with index j = 1, . . . N. The other consideration is that we 
wish to prove surjectivity using only the Jacobi fields ^ sR , because the others are considerably 

more difficult to use in the nonlinear analysis. In fact, using the Jacobi fields R would force 
us to reposition the singularities or change the Delaunay parameters, which for technical reasons 
it is far better to avoid. On the level of linear analysis this is a good choice as well, for with this 
restricted domain the linearization is an isomorphism rather than just surjective. 
We come to the main result of this section. 

Proposition 5 Fix R £ R + , R ^ 1, and /i G ((4 — N)/2, 2) and define the deficiency space 

W = span{^i'+:j = 0,...,N}. (23) 

Then there exists an Sq > such that for all e £ (0,£o]t the operator 

L £ , R : Cj«(£(0, 1) \ {0}) © W > C°^ 2 (B(0, 1) \ {0}) (24) 

is an isomorphism. The inverse map will be denoted L~ X R . In particular, for f £ Cj 3 [ '_ Q 2 (i?(0, 1) \ 
{0}) there exists a unique solution w of (j7^) which has a decomposition 

1 N 

3=0 

The operator 

G £ , R : C°^ 2 (B(0, 1) \ {0}) - C^(B(0, 1) \ {0}) 

and functionals 

K R :C^ 2 (B(0,l)\{0})^R 

satisfy the following properties: 

• G £ _ R is defined for all \i £ ((4 — N)/2, 2) and is bounded independently of e when 1 < /j, < 2. 

• The restriction ofG e .R to the space of functions h(r,9) with eigencomponents hj(r) vanish- 
ing for j ~ 0, . . . , N, is defined and bounded independently of e when —N < fj, < 2. 
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• The functional K® R is defined for all /x > (2 — N)/2 and is bounded independently of e 
when < \l. 

• The functionals K 3 e R , j = 1,...,N, are defined for all ^ > (4 — N)/2 and are bounded 
independently of e when 1 < /i. 

The proof of Proposition [| will be broken into several steps. The solution is constructed by 
separation of variables, restricting the problem to each eigenspace of the Laplacian on the sphere. 
First, we note that although the parameter R is important later, at this stage it is irrelevant and 
may be scaled away by setting y = x/R. In fact, replacing f(x) by R~ 2 f(x/R) and w(x) by 
w(x/R), we see that ( fl9| ) is equivalent to 

f L e w = f in B(0,l/R)\{0} 

\ 25 
[ w = on dB(0,l/R). 

When we change to the variable t = — log \x\, this rescaling has the effect of replacing v £i n(t) = 
v e (t + \ogR) by v £ (t). The K° e R will be denoted simply K§, and so on. In the (t, 9) variables, 
the Holder spaces are converted to the ordinary (translation-invariant) spaces, with the weight 
function e* M . 

Write the eigenfunction decompositions of w and / as 

00 oo 

w(x) = log and f(x) = \x\-^ X>(-logM)&(0). 

3=0 3=0 

If w solves ( poj ) then Wj(t) solves 

MjW = - (— Y~ ) W 3 ~ X 3™3 + ~ ^ L Ve ^3 = U 

Wj (log R) = 0. 

The L CJ - are the restrictions to the eigenspaces of Aj»-i of the operator 

, , d 2 w ,JV-2. 2 . N(N + 2) jX. 
h e w{t, e ) = ^T~ i ^^2~ ) + A s N - lW + — ^ lv " w > 

which is simply problem (H) transformed from M. N \ {0} to the cylinder R x S^ -1 . 
We may assume, after multiplying by a suitable factor, that 

||/||0,o,M-2 = !• 

We first study the operators L e j for j > TV, then for j = and finally for j = 1, . . . , N. 

Step 1. Suppose that — N < fi < 2. Let / be the projection of / onto the sum of the 
eigenspaces of Ag»-i with j > N, 

f(t,e)= £ fjWMQ- 

j>N+l 

We must solve 

( IuMt,6) = f(t,6) in (logi?,+oo) x S' Ar - 1 

S ( 26 ) 

1 u)(log i?,6>) = 0. 
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To do this we first show that for each T > log R there exists a unique solution of 

L e w T (t,6) = f in (logR,T) x S"- 1 
w T {\ogR,6) = 
w T {T,6) = 0. 

This solution can be obtained variationally. Consider the energy function Et 



£ t (w) 



! ' ' * ■ + \v SN - lW \ 2 - N{N : 2) ^ W ^ + J W ) due. 



W 



Since Xj > Xn+i = 2N, we estimate 



£t(w) > 



f J (" ! 

JlogRJS"- 1 V 



: , , ((N + 2f N{N + 2)^\~ 



w z + fw dtd9. 



Since < v s < 1, this energy is convex and proper, and the existence of a unique minimizer for 
Et is immediate. 

Lemma 1 For —N < fi < 2, let 8 — + fj,. Then, for all R > 0, there exist constants e > 
and C > independent ofT, such that for all e G (0, eq], we have 

sup sup e 5t \w T {t,0)\ < C. 
ees"- 1 *e[iog r,t] 



Proof : We prove this by contradiction. Recall that, by assumption, 

sup sup e st \f(t,6)\ < 1. 
e€S JV - 1 te[io g ii,T] 

If the assertion were not true, then there would exist sequences of functions /j, numbers Tj and 
Delaunay parameters e% and a sequence of solutions such that 

lim ( sup sup e st \wTi{t, 0)\) = oo. 

j^oo 0es jv-i log R< t <T t 

We let i)j = w £i . Now choose G (log i?, Tj) such that 

sup e Su \w Ti (ti,0)\ — sup sup e St \w Ti (t, 0)\ = A i: 
ees N ~ 1 ees"- 1 log R<t<Ti 

and define 

Wi(t,6) = A^e^wr^t + U^), 
f i (t,e) = Ar 1 e 5t 'f(t + t i ,e). 

Then, by definition, 

sup sup e st \u>i(t, 6)\ = 1, 

fleS"- 1 log R-ti<t<Ti-ti 

and /j tends to in norm. In addition, 

-^r - 4 w * + A 5«-i^ + — ^ — (* + - fi, 

on [logi? — fj, Ti —ti) x S N ~ 1 . Passing to a subsequence, we may assume that log R — ti converges 
to some number ri G R~U{— oo} and also that Ti — ti converges to T2 G R + U{+oo}. Furthermore, 
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we may also assume that over every compact set of (ti,T2), the sequence Vi(t + ti) converges to 
Uoo(t) and uii converges to w. By construction, w ^ 0. Then Voc is a solution of @ and w satisfies 
the equation 

d 2 w (N-2) 2 „ A . iV(iV + 2) ™i,„ „ 

-jto--- — l w + A S n-iw + — K — '-v£- 2 w = 27 

at z 4 4 

on [ri,T2] x S . Let iJoo = H(vqo, Voo)- There are a few cases to consider, depending on the 
values of t\, ti and H^. To analyze these, we require some growth estimates. 

Lemma 2 For every 77 > 0, iftere exists Eq > suc/i t/iat, /or oZZ j > iV+ 1 and /or a// s G (0, £0], 
every solution of 

^-(^ + A,-^„,")» = ,2S, 

/ (iV + 2) 2 \l/2. . / (W + 2) 2 n \l/2. 

either decays to faster than e ( J '- 1 a£ 00 (respectively, e ( 4 17 - 1 at — 00 J or 
blows up faster than e^ { 4 2t ^ 1 1 at 00 (respectively, e~( ' « 2?? ) X * at — ooj. 

Proof : Because j > iV+ 1, the term of order zero, ^ JV ~ 2 - > + Aj — EMLt3l v ^- 2 j s positive, hence 
L ej - satisfies the maximum principle. Let 

. f{N-2) 2 x iV(iV + 2)\ 1/2 , a . . f{N + 2) 2 ^ 1/2 

*i = P 4 + ^ - " J and 0(f) = P 4 - '? 

We shall show the existence of a solution bounded above by e~^ 2?? )* as t tends to 00 and bounded 
below by e~ /3 ( 2,) ) t as t tends to —00. Replacing t by —t yields a solution with the appropriate 
exponential decay at —00 and blowup at 00. Since these span the space of all solutions, this will 
prove the lemma. 

Since e~ Sjt satisfies h e je~ Sjt < 0, it is elementary that there exists a function yj(t) such that 
yj(0) = 1 and 

Vj{t)>e~ 5it for t<0, <%•(*)< e - *** for t > 0. 

Next, note that 

V& + T *) = yf*T e % yj ® foraU *- (29) 

To see this, simply observe that w = yj(T e /2)y.j(t) — yj(t + T s )yj(—T £ /2) solves (p8|), decays 
exponentially at 00 and takes the value at — T £ /2, hence by the maximum principle must vanish 
identically. We next claim that 

< Vj(t) < e-WiJt for te [ ,T e /2], and Vj (t) > e^ 2 ")* for te[-T e /2,0]. (30) 

Granting this for the moment, then 

v (T /21 e -l3(2 V )Tj2 

Vj(-T s /2) - 



Using this, along with (|29|) and (|30|), the desired estimates for yj(t) for all t are immediate. Hence, 
to prove the lemma, it will suffice to examine the function yj on the interval [— T e /2, T e /2]. 
To prove the claim, we use the fact that 
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To see this, multiply ( ]2q ) by e ljt and integrate from to +00. Integrating by parts and using 
that yj(0) = 1, we get 

Since the right side is nonnegative, we first get that y'-(0) > — 7j- Next, from Proposition ||, we 
can bound v e (t) < 2ee E T 2t for t £ [0,T e /2] and v e (t) < 1 for t > T £ /2. Hence 

JO JT c /2 

Since T e ~ — jv^2 l°S e ' the right hand side of this expression tends to as £ tends to 0, and the 
proof of the claim is complete. 

We finally prove the sharper estimates ( |30| ) on [— T e /2, T £ /2] . From Proposition || again, there 
exist e > 0, t > such that for all e g (0, £ ] and t e [-T e /2 + 1 , T £ /2 - t ] 

(N-2f N(N + 2) (N + 2) 2 2 

1 + x j 1 v e (i) > t] = (3(ri) . (32) 

Using ( |3l| ) we may also assume that 2/^(0) < — (3{rj) for all e £ (0, Eo]- By (|32|), we may use the 
maximum principle on the interval [0, T £ /2 — to] to conclude that 

< y 3 {t) < e^ Wf , < t < T £ /2 - t , and y 3 {t) > e - ^', -T e /2 + t a <t<0. 

Furthermore, y 3 is monotone decreasing, so its maximum on any interval is attained at the left 
endpoint and its minimum at the right endpoint. Hence, for e sufficiently small, 

e -/%)(Tf-to) < e -/3(2r/)^ 

This completes the proof of (|30| ) and hence the lemma. □ 
We now return to our proof of Lemma ^. 

• First consider the case where Hoc ^ 0. Decompose w as 

w(t,9)= «i 

j>N+l 

(N + 2) 2 

If Ti = -co then, taking 2rj < S 2 in Lemma || and using the fact that w is 

bounded as t — > —00 by e~ St with 6 < — K~i we see ^ na ^ ^ or £ sufficiently small, Wj decays 
exponentially as t goes to -co. Similarly, if T2 = +00, then by the same argument uij decays 
exponentially as t — > 00. If t\ > —00, then Wjiji) — and if ti < 00, then Wjfrz) = 0. 

Therefore, we can multiply the equation satisfied by Wj by Wj itself and integrate by parts. 
There are no boundary terms because of the exponential decay and vanishing boundary 
values, so we obtain 

= y (WjH Wj + XjWj Woo u;|)at. 

But since Xj > 2N and < Voo < 1, we conclude that w = 0. 

• Next consider the case where = 0. This case is almost identical to the last one, except 
that the exponential decay is simpler to obtain. Decompose w as before. If t\ = — 00 then 



17 



by Proposition H and (|27j ) Wj grows asymptotically like e ±7j *, where 7? = — 4 2 ^ + Xj, near 

—00. But w is bounded by e~ St near —00, where 5 < — ^- — < 7j, and so uij decays like 

e 7it at — 00. Similarly, if 72 = +00, Wj decays exponentially as t goes to +00, and as before, 
103(71) = and Wjfo) = if either 71 or 72 are finite. 

The same integration by parts as before implies that w = 0. 

In either case, we have produced a contradiction, since by assumption w is non trivial. The proof 
of Lemma |j] is complete. □ 
Finally, we let T —* +00 to get the existence of a unique solution w to ( p6|) which is uniformly 
bounded by Ce" 5 * for t > 0. 

The second and third steps of the proof of Proposition |^ are rather similar to Step 1 above. 
The main difference is that an extra step is needed to ensure that w has vanishing boundary value 
at t = log R. This is where the nonuniformity of as e — > arises. 

Step 2. Let j = and assume that fi > 0. We shall now solve the equation 

f ^s,ow (t) = f in (log i?, 00) 

S 33 
( too (log R) = 0. 

Normalize /o so that ||/o||o,a,/i-2 = 1, and choose an extension /o of /o to K with the property 
that I |/ol |o,«,/i— 2 < 2- For each T > logi?, let wt be the unique solution of 

L e ,owt = fo, w T (0) = 0, w T (0) = 0. 



N - 2 

Lemma 3 For /i > 0, let S — — h /i. Then for all R > 0, there are constants So > and 

C > independent of T , such that for all e € (0, £0] we Ziawe 

sup e 5t |wr(£)l < C. 
te[iogfl,T] 

Proof : First, the two independent solutions of L E) otu = are the functions $° ,=t = r^ - ^g^, 
and neither these, nor any linear combinations of them decay or grow exponentially. In fact, <i>^' + 
is periodic (of period T E ), while 3>°'~ grows linearly. This last fact may be seen by differentiating 
the equality v e (t + T e ) = v £ (t) with respect to e to get 

*°--(t + r e ) + ^^' + (^ + r e ) = $°--(*); 

this gives the linear growth rate since the derivative of T e never vanishes. (Even without knowing 
that T e is never stationary, it is elementary that no solutions of L Ej o can decay at infinity.) 

The remainder of the proof is nearly identical to that of Lemma If the claim were not 
true, then there would exist sequences /o,i, Tj, £j and WT t such that A{ = supjg/^^i e** | iwti | 
tends to infinity. Notice that WTi grows at most linearly as t tends to -co, hence the supremum 
in the last fomula is always achieved. Choosing tj S (logi?, Tj) to maximize e St \wTi(t)\ (so that 
e Sti \wTi(ti)\ = Ai), we rescale the functions and translate the independent variable by ti to obtain 
a solution of 

d 2 w, (iV-2) 2 . N(N + 2) ^ . ~ 
-^2 4 ^ + 4^ u i w i = /o,i, 

in [—00, ^ — tj]. This solution satisfies 

sup e**|tBi(t)| = 1, 
tel-DQ.Ti-tji 
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while /o; tends to zero in norm. 

Passing to a subsequence, we obtain a nontrivial solution w of the equation 

d 2 w (N-2) 2 „ N(N + 2) * „ „ 

over (— oo, T2]. 

• We can immediately see that T2 = 00, for if not then w would be a non trivial solution of 
( j34|) such that w{t2) = w{t2) = 0, which would imply that w = 0. 

• We can also easily rule out the case where H^Vqo, Vqo) = ^ 0. For if w were a nontrivial 
solution of ( |34] ) then it would decay exponentially at 00. But we have already seen that any 
solutions of L £j ow = decays exponentially. 

• The hnal case, where = 0, is also easy to rule out. By Proposition || and equation 
(pi}), to grows asymptotically near 00 like e ±— 2~~ *. Since w is bounded by e~ 5t and because 
<5 > ^j^, we may conclude that w = 0. 

With the proof of Lemma ^ complete, we let T — > 00 and obtain a unique solution to to ( |33| ) 
which satisfies 

sup e 5 *|u)0)| < C. 

t£[log_R,+oo) 

There is no reason why this solution should satisfy the boundary condition at t = log R, and so 
we must add an additional term to correct the boundary data. We define the solution wq by 

w (t) = w(t) - w(logR)(^ + (\ogR))- 1 <S>° £ > + (t). 

This obviously satisfies the equation and the correct boundary conditions, so the proof will be 
complete if we show that i<J>°' + (log i?) is bounded from below by some constant independent of 
e. The following result is deduced easily from Proposition |3[ 

Lemma 4 Given R > 0, R ^ 1, there exists a constant tuq > such that for all e € (0, £0], we 
have 

$°'+(logi?) > em . 



Remark 3 The definition of wo is certainly not unique, for we could have added any appropriate 
combination of and to correct the boundary data. However, is more difficult to 

use in the nonlinear analysis, so it is important that we only use here. Moreover, the use 
of ^g' - would oblige us to change the Delaunay parameters of our approximate solution in order 
to get a solution of the nonlinear problem. 

Step 3. Finally, suppose that 1 < /j,. We wish to solve, for j = 1, . . . , N, the problem 

L £; jtOj(t) = fj in (log R, 00) 

(35) 

Wj(logi?) = 0. 

This is done, mutatis mutandis, following the proofs in Steps 1 and 2. We describe the minor 
changes that need to be made. First, define an extension fj of fj to all of M, and find the unique 
solution wt of 

h e ,jWT=fj, w T {T) = wt(T) = 0. 

The uniform bound of e 5t |tOT| on (— 00, T] is proved as before by contradiction. The fact that 
that the supremum is achieved follows as before from the fact that, for t < \ogR, wt is a linear 
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combination of <j>\' and thus blows up as t tends to — oo at most like a constant times e In 
this step we obtain, by rescaling and translation, a nontrivial solution w of 

d 2 w (N-2) 2 „ , N{N + 2) ^ _ n 

- 4 ' w -{N- l)w + 4 ' v£- 2 w = 0. 

To rule out the existence of a nontrivial solution to this equation, we reason as before. In the 
second case, where ^ 0, we use the two explicit solutions of h Bi jiv = given by ^^(t). 
These solutions do grow at an exponential rate, but by hypothesis, both either blow up or do not 
decay quickly enough, and we conclude as above that this case never occurs. When = we 
note that the asymptotic behaviour of solutions near ±oo is now given by e ~ . 

Now let T — » +oo to obtain a suitably bounded solution of ([35]). Again we must correct the 
boundary data, so we define, for j = 1, . . . , N, the solution Wj by 

w 3 {t) = w(t) - ] ^tD(logii)(^' + (logJ2))- 1 ^+(t). 

Lemma 5 Given R > 0, there exists a constant mi > such that for e 6 (0, eq] and for any 
j = 1, . . . , AT, we have ^^(logi?) > m\£. 

Proof : Recall that 
By Proposition ||, we get 

^ A _1_ / N N — 2 N f _ . W + 2 JV t . 

$l' + (t) = e — - — e~^ 1 + 0(e~ e^*). 

and the result follows. □ 
Again we remark that it will be important in the nonlinear analysis that we have obtained 
a solution with a < J ) ^ + , but no component, for the latter would force the location of the 

singularities of the exact solution to be different from those of the approximate solution. 

Now that we have obtained a solution of the original equation in Proposition || and a (weighted) 
L°° bound for this solution, estimates for the full Holder norm follow by standard scaling argu- 
ments. This ends the proof of the Proposition. □ 



7 A model for the linearization about Ug(R, a, x) 

In the next two sections we finally study the linearization of M about the approximate solu- 
tion Ug(R,a,x) associated to some fixed singular set S = {x\, . . . , X n }, and parameters R = 
(R\, . . • , Rn) an d admissible e — (e%, . . . , e n ), where each Ri > 1. This last condition can always 
be assumed since we may dilate our problem by some factor to n > 0, this will change the set S 
into {kxi, . . . , Kx n } and will change the parameters Ri into As usual, the main point is 

to analyze this operator uniformly as e tends to zero. We shall do this in two steps, first studying 
the somewhat simpler operator 

C SjR w = Aw+ N{N + 2) J2(X(* ~ Xi)u^((x - Xi)/Ri))w (36) 

i=l 

where the displacements Oj have all been set to zero, and then, in the next section, treating the 
true linearization as a perturbation of C s r. Here % is a smooth cutoff function equalling 1 in 
B(0, 1), vanishing outside B(0, 2), and taking values in [0, 1] in the annulus B(0, 2) \-B(0, 1). The 
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point of doing this is that C g r is much simpler to study because the term of order zero is radial 
in each B(xi, 1). 

The main result is that on suitably weighted function spaces, C s ^ is invertible, with inverse 
blowing up as e tends to zero, but in a manner which we can control precisely. We construct the 
inverse for C s r by solving the equation 

C e ^w = f in R"\E (37) 

in three steps: first solve the homogeneous Dirichlet problem for this equation outside the union 
of the balls B(xi, 1); next, solve the homogeneous Dirichlet problem for this equation in each of 
these balls; finally show that the sum of these solutions can be modified to a true solution using 
the Dirichlet to Neumann maps on the boundaries of these balls. 

The main result of this section has a statement parallel to that of Proposition |[ 

Proposition 6 Let e — eq be an admissible set of Delaunay parameters. Suppose that R = 
. . .R n ), n > 3, is a collection of numbers, with each Ri > 1, and satisfying (fH^. Suppose 
also that 

1 < v < 2. (38) 

Define the deficiency space 

Wo = spanixnfi^ 1 ), j = 0, . . . ,iV}. 

Then for e sufficiently small, the operator 

C e -, R : C%(R N \ E) © W — C°!_\_2( kJV \ S ) ( 39 ) 

is an isomorphism. In particular, for each f £ C®'^ 2 -2^M N \ ^) there exists a unique solution of 
fflp which has a decomposition 

N 

"'( 1 " 1 ^ ". 

3=0 



1 

>(*) = GeMM*) + E X(X ~ Xi) & A .(f)-H: + ((x - Xi)/Ri), 

The operator 
and functionals 



G g ,R ■■ C°L a 2 ,_ 2 (M JV \ E) — Cl>«(R N \ E) 



C s .:C!lV 2 (1 W \S)^1, j = Q,...,N, i = l,...,n 
are bounded independently of e. 

We shall denote the right inverse of C s r constructed here by As noted earlier, the proof 

will be done in several steps. 

Step 1 : The exterior problem. Let Q = R N \ U? =1 B(xi, 1). 

Proposition 7 Let f G C u 'f 2 -2^- N \ Then for eq > sufficiently small, there exists a 
unique solution w G C v ' 2 (^) to 

{Us rW — f in Vt 
(40) 
w — on dft, 

where f is the restriction of f to O. This solution satisfies the estimate 

|M|2,Qi,i,,2 < c ll/ l|o,a,z,-2 -2, 

with some constant c > independent of Ei £ (0,£q]. The norms here are taken in fl. 
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Proof : By Proposition |3| the term of order zero in C s r has compact support in f2, and is 
bounded by Ce. We can transform ( ^p|) to a problem of the form Aw' + eV(x,e)w' = f in 
some compact set f2 with smooth boundary by taking the Kelvin transform at some point outside 
of fl. The function V(x, e) is bounded in f2 independently of e. f £ Cl'" 2 -2^) implies that 
/' G C°' a (fl \ {0}). The existence and uniform estimate are now standard. □ 

Step 2 : The exterior Dirichlet to Neumann map. We now introduce the exterior 
Dirichlet to Neumann map. Let * = {ipi(6) , . . . , ip n (0)} G ®" =1 C 2 ' a (dB(x t , 1)) be a set of 
boundary values. It is standard, using Proposition^, that there exists a unique solution C^(51) 
of the homogeneous problem 

rW = in fi 

(41) 

w = ipi on dB(xi,l) 7 Vi = l,...,n, 

which we will denote by w^. The correspondence from ^ to is continuous for each e and R. 
Define S s ^ by 

5 g , fi (*) - (a riW *| ri=1 , . . . , d r „w*\ rn=1 ) G ^ =1 C^ a (dB( Xl , 1)), 
= |a; — Xi\. This is the Dirichlet to Neumann map for the operator C g ^ on Q. 
Lemma 6 The norms of the mappings 

(B? =1 C 2 ' a (dB(x h 1)) 3 ^ — > G C*; 2 a (Q) 

and 

S^h : ©2 =1 C 2 -«(SB( a:i ,l)) — C=iC 1 '"(5i?(x l ,l)), 

are bounded independently of e provided each £, < £o and £o * s sufficiently small. Furthermore, if 
all Si tend to 0, S E # converges to a limiting operator Sq which is the Dirichlet to Neumann map 
for the Laplacian, A, on CI. 

Proof : The boundedness follows from Proposition [?]. Convergence of the S s ^ is a consequence 
of the fact that C E tends to zero uniformly with e. □ 

Step 3 : The interior Dirichlet to Neumann map. We have already proved the analogue 
of Proposition ^| in Proposition ^, so we may pass directly to the corresponding interior Dirichlet 
to Neumann map. As above, the unique solution of the homogeneous problem 

£e.R w = in B(x u l) 

(42) 

w — ipi on dB(xi,l), 
with 

w = w 1 p i G Cl< a (B{ Xi , 1) \ {Xi}) ejLo Span{-^-+((- - x^/R,)} 
defines a continuous map 

* G (B? =1 C 2 > a (dB( Xl , 1)) — > (wfr ,..,%)£ ©r=i K a (B(x l7 1) \ { Xi }) 

©f =0 S P an{i^+((. _ Xi )/Ri)}). 

Now define the map T E ^, by 

r^(*) = (^(^), . . . , Tg R (1>n)) = (d ri w^ | ri=1 , . . . , d rn w^ | rn=1 ) G e? =1 C 1 '"(9S(x i , 1)). 
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Lemma 7 The norms of the mappings 

C 2 ' a (dB( Xi) 1)) 3 — ^ G Cl ,a (B(x l , 1) \ {a:,}) ®f =0 Sfc»n{itf*+((. - Xi )/Ri)} 

Ei 

and 

%,r ■■ ®UC 2 ' a (dB( Xl ,l)) — > e2 =1 C 1 ' a (flfl(x i ,l)) 

are bounded independently of Ei provided all Ei < Eq and Eq is sufficiently small. The mappings 
T s £ converge, as e tends to 0, to a limiting operator Tq which acts diagonally on n-tuples. In 

terms of the eigenf unctions <pj for the Laplacian on dB(xi, 1), the i component T l is determined 
by 

(2-N f(N-2) 2 \ 1/2 \ 
Wi) = ( -j- + [ 4 + *i J j 0, V j > 1, 

T*U^ jV-2 / co S h((iy-2)log^/2) \ / i?f- 2 ^ 



sinh((7V - 2) logiJi/2) ) ^ v \ i?f- 2 - 1 



Proof : We only need demonstrate the convergence of the T E ^ and the specific form of the limit 
Tq. This will follow from the proof of Lemma ^. 

By Lemma ^, for each j > N + 1, there exists a unique solution j/j (t) of the equation 

A ( {N-2f N(N + 2) ^\ 

M {— — + i i £i,Ri )y = 0- («) 

which satisfies yj(0) — 1 and which decays exponentially as t — > +oo. In this equation v Si ^R i = 
v Ez (t + log Ri), Let 7 j = + Aj) 1 / 2 and ^ = ( 7 2 - jv( ^ +2) ) 1 / 2 . As in the proof of that 

Lemma, we show that 

\Vj(t)\ < e-W for all t > 0. 
Integrating the product of ( |43| ) with e~ 7jt from to +oo and integrating by parts, we find that 

dVj , nS , _ /" +0 ° iV(JV + 2) 



^ (0) + lj = J Q ' 4 ' v£* y, (t) e-** dt. (44) 
Since \v Eu R i (t)\ < 1, we obtain 



Moreover, using (44) we see that, for fixed j > N + 1, 

l im ^(o)=- 7 ,. 

Now, the solution of (^) with ipi = cf>j is given by w(x) = \x — Xi\^~ yj(— log \x— Xi\)(J)j(6), and its 
Neumann data is {^^- — ^(0))^>j(#). The relevant solution when Ei = is |x — Xi| 2 ~^ +7i (^(0). 
Hence if ip(9) = J2^^n+i a i < / , j(^)> we see that 



and this expression tends to with e,. 



■7i 
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When j = 0, . . . , N, we can proceed by direct computation. Indeed, when j — 1, . . . , N and 
is the eigenfunction </>.,-, then 



(*^+(logi? 4 )) 1 log \x - Xi\ + l0gRi)4>i(fi) 

is the solution of p2|). An explicit computation, using Proposition ||, shows that 

as expected. Similarly, when j = 0, the explicit solution of (|42j) is given by 

(^ : +(logi? 4 )) _1 ^: + (-log|x-^|+logi? 4 ). 
Using Proposition || again, we get 

,. ~i u , N ~ 2 (i ■ cosh((iV-2)logfl i /2) 
= — ^ + S mh((iV-2)log^/2) 



□ 



Step 4: Invertibility of the difference of Dirichlet to Neumann maps. In order to 
glue the interior and exterior solutions together, we must add correction terms to these solutions. 
These correction terms are solutions of the homogeneous problem, and are chosen so that the 
Cauchy data from the inside and outside match up. To find these, we must show that the 
difference S g ^ — T E ^ is invertible when all £j are sufficiently small. 

Proposition 8 There exists £o > such that for e < Eq, 

S S ,R ~ %M ■■ ®tiC 2 ' a (8B( Xi , 1)) — > ®U^ a {dB{ Xl , 1)) 

is invertible. The norm of its inverse is bounded by some constant C > independent of e. 

Proof : By the L 2 operator norm convergence of S g ^ — T s ^ to Sq — T , it suffices to show that 
S — T is invertible. Now <S — % is a self-adjoint first order pseudodifferential operator. Both <S 
and T are elliptic, with principal symbols |£| and — 1£|, respectively, hence the difference is also 
elliptic and semibounded. This means that So — To has discrete spectrum, and thus we need only 
prove that it is injective. The invertibility in Holder spaces then follows by standard regularity 
theory. 

We argue by contradiction. Assume that So — To is not injective. Then there exists some 
*o G ®? =1 C°°(dB(xi, 1)) for which (5 -7o)(*o) = 0. We may extend the Dirichlet data * to a 
harmonic function on the exterior region, decaying at infinity, and also one on the interior region. 
So is the Dirichlet to Neumann map for the ordinary Laplacian on the exterior region, while To 
is the Dirichlet to Neumann map for the Laplacian on the union of balls only for the spherical 
eigencomponents with index j > 1. Thus \P extends to harmonic functions w' and w" on the 
interior and exterior regions, respectively. These functions have the same Dirichlet data, namely 
an d by assumption on this Dirichlet data, all eigencomponents with j > 1 of their Neumann 
data agree as well. To find a harmonic function with the j • = eigencomponent of the Neumann 
data also matching at these spheres, it suffices to add on appropriate multiples of the functions 
\x — Xi\ 2 ~ N — 1 on each B{xi, 1). We obtain in this way a harmonic function w on M. N \ E which 
decays at infinity, and with at worst a radial singular term at each a;,-. The only possibility is 
that this function has the form 

n 

w(x) = J^Pilx - x l \ 2 ~ N . 
»=i 
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We now find the radial part, Wq(x), of the harmonic function w(x) in each ball B(x^, 1), 
and use this to compute Tq^q directly. Comparing this with the previous answer will lead to a 
contradiction. In B(xi , 1), this radial part is a sum of the singular part and a constant: 



m UJN - 1 J*"- 1 



-1 Js" 



Here ojn-i is the volume of S N 1 . The terms in the sum may be computed by noticing that for 
any fixed vector v, the function 



1 



\x + v\ J -"dO 

is harmonic, regular and radial in B(0, \v\), hence constant. Thus, provided |u| > 1, its value 
when \x\ = 1 is the same as its value at x = 0, and is just ul 2- ^. We conclude that 



w l ° (x) = pi a \x - x io | 2 N + ^2 Pi\ x i - x 



2-N 



The normal derivative of this function at B(xi , 1) is just (2 — N)pi , but on the other hand, it 
must agree with Tq° (uJq° ) . This leads to the equality 



2-N 



R N ~ 2 I — , 
(2 - N)p io = (2 - N) 1" 2 - p io + Pi\xi ~ Xi 1 

or simply 

ft +<- 2 E^i^-^oi 2 - iV = o. 

Now recall that the parameters Ri satisfy the relation (|l2]), 

^2(Ri)^(R io )^qi\x io ~x,\ 2 - N = q io . 

Pi / 2-JV 

A small calculation shows that the collection of numbers pi = — (Ri) 2 give a solution of the 
system 

Y(Ri)^{Ri o ) i£ ^'qiqi o \xi o -Xi\ 2 - N (p i +pi o )=0, i = 1,..., n. (45) 



JV-2 . _ . N— 2 



Let Sjj = (R4) 2 (Rj) 2 qiQjlxj — Xi\ ■ Multiply ( |45| ) by pi and sum over iq to obtain 

n 

S hio(Pi +Pio))Pia = Y s i,i (pi + Piof = °- 
io=l «#»o i<«o 

Since all the Sj,^ are positive, we get that p, = — Pi . But since n > 3, this forces the vanishing 
of all pi. This is a contradiction, and the proof is finished. □. 

Step 5: The correction term. Let / G C®'" 2 -20"^ \ ^0- By Proposition [| and Proposi- 
tion ^ there exist functions w ex t and Wi n ^i, i — 1, . . . , n, such that 

J £e,RWext = f m fi 

i w eKt = on <9fi, 
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and, for i = 1, . . . , n, 

( C St RW in t,i = f in B(xi,l) 
\ w int} i = on dB(xi, 1) 
From Proposition B, there is a decomposition inside each B(xi, 1) 



w int ,i(x) = G euRt (fi)(x) + f^Ki uRi {f i ){x)j^+{^i), 
where /, = / Isfe,!)- I n addition, for some c > 



JV 

\\G El MM)\ka,u + Yl \ K LnM*)\ ^ c||/|| ,«,„-2,-2 
J'=0 

and 

||«?ea;t||2,o,v,2 < c\ |/| \o,a,v-2,-2- 

We now seek a function Wker in K w \ {Uf =1 dB(xi, 1) U £} so that 



w(x) 



W ex t{x) + W k er(x) in SI 

Wint,i(x) + w ker (x) in i = l,...,n, 



will satisfy C^ R w = / on R N \ S. 

For this to be true we need to choose Wk er to be continuous across each dB(xi, 1), to satisfy 
C s RWker — away from these boundaries (and from S), and such that the jump of d ri w^ er across 
dB{x il 1) equals d ri Wi nt .i — d ri w ext - This is equivalent to finding a solution '5 of 

(S gt& - T g)R )($) = {(d ri w mt ,i - d ri w ext ), (d rk w int ,k ~ d rk w ext )} G ®^ 1 C 1 ^ a (dB(x l , 1)). 

This is possible by Proposition ^ this same result gives the uniform bound on the norm of ty. 
All estimates for the norm of w follow from Lemma ^| and Lemma ^. □ 



8 The true linearization 

We now undertake the main task of analyzing the linearization about Ug(R, a, x). Recall that our 
ultimate goal is to find a solution of the nonlinear equation 

N(N-2) - n+2 
j\f(w) = A(u e (R, a, •) + to) + — ^ L (u g (R, a, •) + w)— = 0, (46) 

for some choice of parameters a,R (recall that s is regarded as fixed), and for some w which is 
dominated by Ug(R, a, x) near £ . Let 

A sA - a = A + N{N + 2) (u e -(R, a, -))^ (47) 

denote the linearization of this equation around w — 0. As is usual for a Jacobi operator, one can 
obtain solutions of A g ^ a w — as derivatives of families of solutions of (^6|). Of course, we do not 
know any families of solutions yet, but we can obtain approximate solutions of A g p, a w — from 
families of approximate solutions of ([Iff), and we know many explicit families of approximate 
solutions. In particular, the functions 

9 g',R,a,i^ = dRiUg(R,a,x), and _ . (x) = d a jU g {R, a, x), j = l,...,N, i=l,...,n, 



26 



defined by differentiating Ug(R,a,x) with respect to Ri and a\, respectiveiy, have the property 
that 

vanishes near each Xi . This may be seen by differentiating (^) with respect to either Ri or a\ . 

In this section we shall consider the problem A g j^ a w = /, for / S C°l? 2 _ 2 (K W \ £). The main 
linear result in this paper is that we can find a solution with optimal bounds in the space 

M{a) = C%(R N \ S) © Span {^;+ a J, (48) 
where the norm of an element w(x) = v(x) + Y2i a i Si SjLi a i % a i is given by 

IklU(a) = P v \\v\\2,a,u,2+S^2 \ S i\ + £ P l^l" 

As before, the deficiency space spanned by the a i is necessary to obtain surjectivity of the 
linearized operator. The rather byzantine choice of norms on M.(a) is necessitated by the rather 
intricate estimates below, and the need for uniformity as e tends to zero. The following proposition 
states the expected bijectivity of p, a , but it also gives a rather sharper estimate on the solution 
in terms of the linear functionals KL 3 - R . introduced in Proposition^. The values K, J - Ri (f) are 

closely related to the coefficients Si and ol\ for the solution w £ M(a) of A s ^ a w = /, but they 
correspond to coefficients of the simpler, radial solutions . It may seem rather unnatural to 
continue to use these radial solutions here, but once again, these sharper estimates seem necessary 
later. 

For notational convenience, we shall often denote these approximate solutions by 

( ''" /,. „. , {x)=ih (e, R,a,x), * |± _ . (x) = (e, R, a, x), (49) 

and we shall identify w{x) = v(x) + SiPi(e, R, a, x) + a^f (e, R, a, x) with (v, S, a). 

Proposition 9 Let e be an admissible set of Delaunay parameters, and suppose that Ri > 1 and 
ai G l w are a collection of numbers and vectors satisfying and (j7j|j. Then 

is surjective, provided all are less than some sufficiently small number Eq. Furthermore, for 
f 6 -i^ N \ there exists w = (v, S,a), corresponding to 

N n N 

w(x) = v(x) + ^^Si^i(e,R,a,x) + ^ ^ oc^jj (e, R, a, x), 

i—l i=l j=X 

solution of the equation A g R a w — f , which satisfies 

p"IMka,,,2 <c(p"|l/ll0,a,,-2,-2+SUp|^ ifiii (/)|+pSU P |4 fi< (/)|). 

e\S\ < c{p v+1 ||/||o, ai „_ 2 ,_2 + sup \K? Ri (f) | + p 2 sup \K? g Ri (f)\), 

» i,j ' ' 

and 

ep\a\ < C ^ +1 ||/||o,a,,-2,- 2 + psup |/C° & .{f)\ + psup \K? g R ,.(/)|). 
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Proof : Define the space 

M = C 2 J(R N \ E) © Span{xOz - Xi)£i,x(* - x^}, 

where now \ is t ne cutoff equalling one in each B(xi, 1) and vanishing outside the union of the 
B(xi, 2) and 

fli(ei,Ri,x) = d Ri u £i (Ri,0,x - Xi), 7- (e t ,Ri,x) = d aJ u £z (Ri,0,x - Xi). 
An element w = u + 5*ix(x — £Ci)/ii + o^xK^ ~~ ^OtI nas the norm 

\\w\\m = \\(v,S,a)\\M = \\v\\2,a,v,2 + £*y] \Sj\ + £^|a||. 

Proposition ^| may then be rephrased as the assertion that 



-e,R 



M — Cl q 2 ,_ 2 (K jv \S) (50) 



is an isomorphism, and that the inverse is bounded independently of e. 

We shall prove this proposition by reducing it to Proposition ^ by a sequence of perturbations. 



Step 1. We first claim that 

A e - b = A + N{ ' ■ ' (u E -(R, 0,x))^:M- -Ck ^ \ ^ W 



is also an isomorphism, with inverse bounded independently of e. To see this, we estimate the 
norm of the difference A = A g ^ — C s ^. First, the difference 

\(u e (R,0,x))T^ ~^2x(x - Xi)(u £i (Ri,0,x - Xi))-^\ 



1 . 

vanishes in each |a; — as« | < Pi and is estimated by cs N ~ 2 dist(ic, E) when dist(x, E) > 1. We get 
the estimate 

eT^(\x-Xi\ N - 6 p+\x- Xi \ N - 5 ) 

for this difference in the annulus pi < \x — Xi\ < 1 from Corollary^. 
Finally, from Proposition ^ we derive the bounds 

\fk(si,R%,x)\ < Ce\x - Xi\ 2 ~ N and pyf (e*, R4, x)\ < Ce\x - Xi\, (52) 

which hold in particular for all x satisfying dist(a;, Xi) 6 [pi, 2]. Thus, for w £ M, we get 

n 

— 4 ^ — ^ 4 

\\(Ag,R,0-£e,R) w \\o,a,u-2,-2 = | \(u S (R, 0, X)"^ - ^ x(x ~ Xi)U ei (R4, 0, X-Xi)^)w\ \o, a ,u-2,-2 

i=l 

< cp N+2 {\\v\\ 2>a , v>2 + ep-'- 1 \S\ + s(l + p N - 2 - v )\a\). (53) 

By (|l0|), we can estimate the right side by an arbitrarily small multiple r\ of ||w||jvi, provided e 
is sufficiently small, i.e. \Aw\ |o, a ,v-2,—2 < ^IIHI-M- Letting £7^ denote the inverse of the map 
©, then 

Kko - + - + 1 (M) 

is a right inverse of the map (|5l|). Since £g ^ is invertible, and since A is small in operator 
norm, A s ^ is also invertible, hence ( p4[ ) is the unique inverse. Since A and £7^ are bounded 
independently of e, so is A7^ Q . 

We now obtain sharper estimates for the solution of A s ^ w = f. 



28 



Lemma 8 Suppose that f is any function in C u '" 2 - 2 > an< ^ ^ w ~ ^e^RO? = ('-'''-'> a ) an ^ 
w a = C~jif = (v , S Q ,a ). Then 

IMIw.3 < C(\\v \\2,a,u,2 + ^""elSbl + p N+2 (l + /^"^laol) , 
e\S\ < C(p N+2 \\v \\ 2 , a , v , 2 + e\S \ + p N+2 (l + p N ^~ v )e\ao\), 

and 

e\a\ < C(p w+2 ||« D || W ,2 + P N+1 - u s\So\ +e\a \). 



Remark 4 We point out explicitly again that the coefficients So and ao here are the values of 
the functionals ICi Ri (f) introduced in Proposition^, corresponding to the radial model operator 
Cg R . Specifically, 

e l (S ) l = /C°^(/), £<(<*,)? = KP gR .{f). 

Although it might seem more natural to view the coefficients S and a coming from the solution 
of A gR0 W — f (and later, K sR5 w — f) in terms of them, it turns out to be simpler in the long 
run to view these as primary. 

Proof: Rewrite ( p4| ) as 

Kko = C t^ + ^RrHl + AC- R - AC , r = CZ^I -(1 + AC J, i ',!£> 

Applying this to / yields 

w = wo - C-\(I + AC-^Awo. 

The estimates of the lemma follow by considering the different components of the two sides of 
this equation, and using the boundedness of £7^(7 + ACZ 1 ^ 1 , and (|53|). □ 
Using these estimates, and the fact that | |w o 1 1 A4 < C||/||o,a,i/-2,-2, we conclude 

1^112^,^,2 < c(||/||o, ct ^-2.-2+p JV+1 -"sup|/CO^ .(/)|+p JV + 2 (l+ / 7 JV - 2 --)s U p SUp |/C^.(/)|), 

i i j=l,...,JV ' ' 

e\S\ < c(p N+2 \\f \\ , a ,„-2,-2 + sup |£° R .(/)| + p N + 2 (l + p N - 2 -) sup sup \K* S R 4 (/)|) 

i i j=l,...,N ' ' 

and 

e\a\ < c(p N+2 ||/|| 0> a,,-2,-2 + p N+1 -» sup |/C° R .(/)| + sup sup \K* R 

i ' ' i j = l,...,N ' ' 

Step 2. We now perturb further. We claim that A. s R is invertible from A4(0) to -2j 
with inverse uniformly bounded in e. In fact, because of the uniform invertibility of (51), it 
suffices to show that the inclusion 

M — ► M(0) 

is uniformly bounded. Equivalently, we must show that 

IMI-M(O) < C||w||a4 

for some C independent of e. To calculate the norm of w in A4o, we must write this function 
in the form w — v + J2 Sifii(e, R, 0, •) + J2 a ilf & 0, •). Since w = v + x(J2 SiPi{et, Ri, •) + 
E^Tife.-Ri.O) we get 



v = v 



+ s i(xi*i(ei,Ri, •) - m(e, R, «, •)) + 4{Xl j M,Ri, - R, a, ■))• 
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Now observe that 

Mi "Mi = (1 - Xi)dRi{Uei{Ri,0,X - Xi) - Wg), *j{ - jf = (1 - Xi)^ («e< (-Rj, 0, X - Xi) - Wg) 



in each B(xi, 1). From (15) 



I IxMi - Mi I ka,^2 < "e, and \\xif - if \U,a,u,2 < Cp 1 "e. (55) 

Hence 

||X)^(XMi-Mi)lka,v,2 + ||^c^'(x^-7i')lk Q ^2 < Cp-"(e|5|+ep|a|), 

and thus 

IMLm(o) = P"||^l|2,a,i/,2 +ep|a| < C(||u||2,a,i/,2 +ep\a\) = C||u;||a4, 

as desired. 

We also derive the precise estimates for the inverse of 

K-,R,o ■■ M(0) — C°'\_ 2 (R N \ E) (56) 

from the ones stated at the end of Step 1. In fact, replace v by v there and use ( j55"| ) and the 
sharp estimates for e|5| and £/c|a| to get 

P"\\v\\2, a ,v,2 <c(p tf ||/||o,a jV -2,-2 + SUp|/C° fii (/)|+/>SU P Sup |/C^.(/)|), 

i i j=l,...,JV ' ' 

< c(p iV + 2 ||/|| , ct , I .- 2 ,- 2 + sxap |/CJ ^ ^(Z)! + y0 JV+2 (l + p^- 2 --) sup sup |AC^.(/)|), 

i i j = l,...,N ' ' 

and 

e P \a\ <c(/ r+3 ||/|| (W ,-2,-2 + p N+2 ~ v sup .(/)|) +psup sup 

i i j=l,...,N ' ' 



Step 3. We perform a final perturbation of (56) to show that 



A ./,..„• : M(a) — C^VaO^ \ E) (57) 

is an isomorphism, with inverse uniformly bounded as e tends to zero. In order to apply the same 
sort of perturbation argument, we would like A s ^ a and A g ^ to be acting on the same space of 
functions. To this end, define 

M = {(v, S,a) e Cl:°{R N \ Z) ® {R+) n ® (R N ) n }, 

with norm given by 

\\( v ,S,a)\\j^ = p v \\v\\2, a , v ,2 +s\S\ + pe\a\. 
Then M. is equivalent to any of the spaces M.{a) via the map 

L a : (v,S,a) — ► v(x) + '^2s i p i (e,R,a,x) + orjljie, R, a, x). 

We denote 

Knra = KRfi °><a:M^ C°^ 2i _ 2 (R N \ E). 
The uniformly bounded inverse A7^ is a first approximation to the inverse of A g j=. a . In fact, 

s,R a A-e,-R,o = -f + ^c fi o( ~ Aj fl o). 
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We shall show that the final term A_ ^ (A g R a ~ A £ R o) nas small norm as a map on A4. Because 
of the uniform boundedness of AT^ n , it should suffice to show that the difference 

K-,R, a ~ K-,R,o ■■ M — C°^_ 2 (R N \ E) 

has small norm. This is almost true, although at one stage of the argument we require slightly 
more information about the inverse. However, let us concentrate on showing that most terms in 
this difference are small. 

Let w = (v, S,a) £ M, and for convenience write 

U a (x) = ^2 SiHi(e, R, a, x) + ^ a Wi & 5 > x )- 

Then 

(K-,R,a- hn,o)( w ) = («g(fl,o,i)^s -u £ (R,0,x)^)v + (Ag^Mx) - A^ Rfi U (x)). (58) 

We estimate the two terms on the right in turn. 

The first is simpler. Going back to the definition of Ug(R, a, x), we expand in B(xi, pi) to get 

(u e (R,a,x))->£* - (u e (R,0,x))^ = 

G-N 

2a, ■ (x - Xi)\x - Xi\~ 2 Vsr 2 (-log|x - x;| + log-R,)(?; £i (-log |x - x,| + logi^) 

2 t; £l (-log|x-x l |+logi? l )) + 0(l)w e 1 T , (-log|x-a; J |+logi? l )- (59) 



N - 2 



Here we have used the fact that for all T > 0, there exists some ct > 0, independent of e < Eq, 
such that for all t S [-T, T] 

v £ (t + t ) + |« e (t + *o)| + |S e (* + *o)| < c T v e {t). 

We consider this difference in three separate regions about each Xi. In the innermost one, 
where dist(x,Xi) < = e N ~ 2 , since v £ and v £ are bounded we get 

\(u e -(R,a,x))^ - (u g (R,0,x))^\ < C\x - x^ 1 . 



For the middle region, where dist(x, xi) € [s N ~ 2 , pi, we further expand (59), using Proposition |3j 
to get 

\{ug(R,a,x))T^ - (u g (R,0,x))Tr=3\ < C{e^\x - x t \~ 2 +e^\x - x^ 5 ). (60) 

A similar estimate also holds if dist(x, xi) S \fii, 2pi\. Finally, when dist(x, xi) > 2pi for all i, then 

(u s (R,a,x))^ - (%(£,(), x))^5 = o, 

since u s does not depend on a if dist(x,Xi) > 2pi. The estimates for the full Holder norms are 
similar. 

Putting these estimates together, we get 

\\(u s (R,a,x)-^ -u g (R,0,x)^)v\ \o,a,v-2 -2 < ce^\\v\\ 2 , a ^,2- ( 61 ) 
The second term on the right in (m) is more complicated to estimate. First write 

A ,; „(' .,) - A= 5 (U ) = A(C/ a - U ) 
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iV(iV + 2) 



4 

Recall that C/ a satisfies 



({u g (R,a,x))^U u - (u s (R,0,x))^Uo) . (62) 



Au _ + N ( N + 2 ) ( Ug (R^,x))^U- a = 

when dist(x, Xi) < pi. Hence in this ball, (|6^) vanishes identically. It also vanishes when 
dist(x, Xj) > 2pi, because Ug(R, a, x) does not depend on a and so the difference cancels. 

We come at last to the transition annulus, where dist(x, x{) £ [pi,2pi\. Recalling that Ug, is 
the differential of ug(R, a, x) and that ug(R, a, x) — Xi(x — Xi)u £i (Ri, a^, x — Xj) + (1 — Xi)^e{x), 
we write 

U- a = XiK + (1 - Xi)W g 
where [/* and Wg are the differentials of u £i and wg, respectively. In this annulus we can estimate 

\U^\<C(p 2 - N e\Si\ + £ P \a\), and \W g \ < Cep 2 - N \Si\. (63) 

Now write ( |62| ) as a sum of three terms, I + II + III, where 

/ = (A Xi )(E£ - U*) + 2V Xl • V(C/; - 

// = MK, -U*)+ X i N{N A + 2) a, x) ^ - Ug(R, 0,x)*±*U$, 

and 

III=(l-Xi) N{N 4 + ^ W g (u g (R,a,x)^ -%(i?,0,x)^). 

In the last of these, we have used that Wg does not depend on a. 
This last term is the simplest to estimate. In fact, by Corollary 

\u £ (R,a,x)T^ - u e (R, 0, x)f^\ < e™^p N ~ 5 . 

Hence 

||//I||o J a,v-2,-2 < Cf?- V ef?- N e*t*p fr -*\S\ < C{p N+x - v )e\S\, (64) 

using @. 

To handle the middle term, recall that 

, N(N + 2) , , 4 , 

^ + —^-4 -^(iZi.Oi.S-SEi)^^ =0. 

Using this, we rewrite // as 

N(N + 2), - _ - • \ 

X—^ {(u s (R, a, x) -u ei (Ri,Oi, x) )Ul. - (u E {R, 0, x) "-2 -u e , (R i ,0,x-x l ) )U^j . 



Since 
we get 



\ug(R,a,x)"-z -u £z (R l ,a u x-x i )N i -2\ <Ce^p N 5 , 
|//||o,a,,-2,-2 < C P 2 - v {p 2 - N e\S\ + ep\a\)e^p N - 5 < V (e\S\ + ep\a\), (65) 



where rj tends to zero with e. 

The first term, /, is the most difficult to estimate suitably. First observe that 

U a z —Uq = SidB,i{Ue,{Ri, a i, x ~ x i) ~ U Ez (Ri,0,X - Xi)) 

+ / ] o4.d n i (u ei (Rj, aj, x - x^ - w £i (i? 4 ,0,x - x^). 
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From the proof of Proposition |j we see that 

N — 2 2 ~ N 

u £i (R l ,a i ,x - Xi) - u ei (Ri,0,x- x{) = e, — - — R. 2 a r ■ {x - x % ) + 0(e\x ~ x t \ 2 ). (66) 

Noting that the derivative with respect to a of the leading term of this expression is independent 
of a, hence cancels in the difference, we get 

\Ui.-Ui\<Cp{e\S\+ep\a\), and | V(W a - W )\ < C(e\S\ + ep\a\). 

Thus we see that both terms in / are of size + £/9|a|), and in particular, they are not 

small relative to the norm of Ai. Fortunately, the term we actually want to estimate, Ar^ (J), 
is still small because I has support in the union of annuli B(xi, 2pi) \ B(xi, pi), which has small 
volume. The following result will be proved in the next section. 

Lemma 9 Suppose v € (1, 2). Let h be the term I . Expand it in terms of the eigenf unctions on 

oo 

h( Xi +rd) =£^(r)&(0), 
j=o 

and let 

N oo 

/i' = ^(r)^(0) and h" = £ fy(r)&(0). 

3=1 ]=N+l 

Then the solution w — C~ R h — (v,S,a) = (u, KP(h), )C'(h)) satisfies 

\\v \\2,a,is,2 < C(p 2_,y SUp|/lo| +p 2_I/ SUp|/l'| + | \h"\ \o, a ,u-2) ■ 

In addition, for any p S (—-/V, 2), we have 

e\Si\ = \ICl R l (h)\ < C( P 2 sup \h (r)\ + P N+1 sup \h'(r)\ + \\h"\\ , a ^- 2 ), 
and, for j = 1,...,N, 

e\4\ = \K.i R>i (h)\ < C(p 2 sup \h a (r)\ + p sup \h'(r)\ + \\h"\\ , a ^-2,-2). 
If i' ^ i, then for j = 1, . . . , N, 

\4\ ee \ICl R i ,(h)\ <C{ P 2 sup \ho(r)\+ p N+1 sup \h'(r)\ + ||fc"||o,«,„-2,- 3 ). 



The point of this result, of course, is that because h is supported in a region with small volume, 
one can obtain better estimates for C-~^ R h (and then also for AT^ Q /i) than the obvious ones. We 
apply this as follows. Using (66), we see that the terms in its eigenfunction expansion satisfy 



N < e{\3\ + \a\), \h'\ < Cp-\e\S\+ep\a\), 



and 



||/i"||o,a, M -2,-2 <Cep 2 -»(\S\ + \a\), 
for any p G (-N, 2). We also see that 

p v \\h\\o, a , u -2,-2 < C(p 2 sup\h a \+p 2 sup | ft' | + p u \\h" \\o, a ,„ -2, -2) < Cp(e|S|+ep|a|). 
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Now apply the Lemma to find a solution £_ = w = (v, S, a). It satisfies 

p u \\v\\ 2 ,a,u,2 < C(p 2 (e(\S\ + \d\) + p 2 (p-\s\S\ + ep\d\) + f(e?-"(\S\ + \d\)), 
s\S\ < C(p 2 e(\S\ + \a\) + p N (e\S\ + ep\d\) + p 2 ^s(\S\ + \a\)) , 

and 

e P \a\ < C(p 3 e(\S\ + \d\) + (?{p-\e\S\ + ep\d\)) + ep^{\S\ + \d\)). 

In the second and third expressions, we have estimated h" in the C '" 2 _ 2 norm, for any p € 
(-N, 2). By taking p = 0, for example, we see that we can bound all three of these expressions 

by 

Cp(e\S\+ep\a\). 

The solution w = AZ ^Ji that we are actually interested in can be estimated in terms of these 
quantities, using Lemma |§|. Combined with the preceding, we see that 

\\w\\m < p{s\S\ +ep\a\). 

We combine this estimate with the estimates ( |6l| ) , ( |64| ) and ([35]) , and use the uniform bounded- 
ness of Q for these terms to finally conclude that the operator norm of A~^ (A g ^ a — A g p . ) 

on M. can be made as small as desired by choosing e small. This will show that 

e,R,0 ^S,R,a = I + A- \ (Ae,R,a ~ ^e,R,o) 

is invertible, hence 

G = (I + A-^ (A sAa - A^o))^-^ 
is a left inverse for A f ^ a . Clearly it is uniformly bounded in e. Direct computation shows that 

KR,a G = ( A s,i?.,o + (A £ -,fl, a - h,R,o))G = Ag Ao (J + A-^ (A e -_R )a - L Rfi ))G = I 

as well. Thus A g R a has a two-sided inverse, which is therefore unique, and finally, A s R a does 
as well. 

Step 4. We conclude the proof of the theorem by showing that the sharper estimates for the 
solution of A g R a w = f are valid. We proceed as in Lemma ^| Thus, let wq = (vq, So, do) = 

KjiflW and w = ( w ' ^' = K.R.J- Then writin g B = A e,R,a ~ ^g,R,o, we derive 

A = v+ A /,•.••''•' = (' + A .;,.,,'<' - A /,•.••'' ' + A S,o s r lA S,o = 

A .L V;,,,/^ • A .«,, /?i ' A J,'::' 

Apply this to / to get 

W = w - Agji fi Bwi, (67) 
where we have set wi = (I + Al^ Q B)~ 1 wo- Clearly we have 

IKII^<C||«;o|Im- ( 68 ) 

We divide w± into a sum v% + Wi, where we have just lumped the second and third terms in the 
usual decomposition into uix. Thus if we let 

~Kji,Q Bvi = ^> 3> and ~ KjifiB™ 1 = ^' ^' 

then we have 

(w, S, d) = (v , S , do) + (v, S, d) + (v, S, a). 
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We estimate these new terms in turn. First, from (|6l|), 

i N + 2 N + 2 

\\A^ Ro B Vl \\ M = p u \\v\\ 2 ,a,u,2 + e\S\ + e P \a\ < Cp-^-\\v 1 \\2, a ,ufi<Cp-^- v \\w 1 \\j Cl . 



To estimate this, use (68) and the sharp estimates for wq from the end of Step 2 to conclude that 



this term is bounded by 

+ sup \K° gtA .(f)\ + psupsup \KP g Ri (f)\). 
I I j ' ' 

For the other term, use Lemma || to estimate ||AT^ ^iSiH^ first in terms of wi and then w . 
Finally, substitute the sharp estimates from the end of Step 2 for wq in terms of /. We omit the 
details, which are straightforward. In the end, we get altogether 

P"\\v\\2, a ,v,2 <c(^||/l|0,a, V -2,-2+SUp|^ 5i (/)|+pSU P SU P \K? gRi (f)\), 

i i j = l,...,N ' ' 

e\S\ < C (^ +1 ||/||o,a,^2,-2 + sup|/C°^(/)|) +(0 2 sup sup 

i i j=l,...,N ' ' 

and 

ep\a\ < C (^ +1 ||/||o, Q ,,-2,-2 + psup|/C^ 2 (/)|)+psup sup \IC gRi (f)\). 

i i j=l,...,N ' ' 

The proof is complete. □ 

9 Estimates on the error term 

In this section we shall derive estimates on the size of the error term in the approximate solution. 
Lemma 10 Suppose that the parameters Ri > 1 and ai 6 M. N satisfy the relations ( fZl| ) and (W 

W+2 

— ; — -u^~ 2 ; then set w = CZ s 

4 e j e.R 



of Proposition 0. Let u g — u s (R, a, ■), and define £ = Au g + N 4 2 ; then set w = C-^Q 

Then, for v £ (1,2), and for some constant C independent of e, 

\\w\\ u <Ce P 2 . 



Proof : Write £- sR C = w = (v,S,a) as usual. The estimate for v is rather straightforward; 
however, the estimates for S and a are more delicate, since they rely on the particular structure 
of C) and so we need to follow the construction of CZ R in Proposition ^| step by step. 
Step 1. From the definition of rig, we see that in 1^ \ U? =1 B(x l , 2 Pi ) 

C(as) = — ^-j '-ur 2 (x), 



and so 



|C(a;)| < Ce^dist(a;,E)" 



In particular, in Q = l w \Uf_ 1 B(xi, 1), |C| < Cs N - 2 . From Proposition 0, there exists a solution 
of 

C g p,w ext = C in fi 



■><> ( 69 ) 
w ea;t = on os 2. 



which satisfies 



JV + 2 
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Step 2. Our task now is to estimate £ in each B(xi, 1). In this step, we consider its eigen- 
f unction expansion, 

+ OC 

C(*) = E (0), 

3=0 

where r — \x — Xi\, and estimate the terms in three separate regions. 

N + 2 

From the definition of u s , we see that inside each B(xi, pi), £ = Au s H 4 ~ Ug~ 2 — 0, and 

so all Q = 0. 

Next, C = N( - N -^ uf^ in B{x i: 1) \ B(x t ,2 Pl ). Thus 

AT 



3 = 1 



7V-2 f ° 



and 

I E CjM^-WI^^^^ 2 - (70) 

j=JV+l 

Finally, in B(xi, 2pi) \ B(xi, pi), we have 

C(x) = 2(d ri Xi)(dr t u ei (R i , at, x - a?i) - d n w s (R,x)) + (Axi)(u El (Ri,a,i,x - x t ) - w s (R,x)) 
_7V(iV — fy^ Xi ( u ^ R .^ auX _ Xi ))w^ _ ( Xl u £i (R l ,a l ,x - x,) + (1 - x%)w g {R, x))%%) . 
From this we get the estimates 

ElCiWI + l E Cj(r)^WI<C £ . (71) 

i=o J>iV+l 
Step 3. By Proposition ||, we get, for i = 1, . . . , n, solutions of the problem 

C gR w int ,t = C in B(x. L ,l) 
w in t,i = on dB(x il l). 



These satisfy 

N 

Wint 



N 1 

Ge.MCm^i)) + E^ A (C|B(, I ,i))-*i' + ((2: - Xi)/Ri 



3=0 

We now estimate the different terms in this formula. 

First consider the solution Wq in the j = eigencomponent. We can represent w l explicitly 
in terms of Q using the variation of constants formula. Write 



&(r) = r^-- 5 ^(-log(r/i2 i )) = r^-*°;+(-]og(r/iZi)). 



2-W 

"at 

Using Proposition ||, and the fact that R t > 1 + c , we see that, up to a constant, ip e ' + R. ~ sr 2 ~ N 
for p < r < 1. In particular, this function never vanishes. This means that we can write 

w l Q (x) = K° iMi (C lB{XtA) (x))^^((x - Xi)/Ri) + wl{x - Xi ). 

where 
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and the first term in the decomposition of Wq is intended to correct the boundary value. From 
the bounds above on Q, and using the approximation above for tp, we find that 

K| < Cep 2 for all pi < \x - x t \ < 1. (72) 

Therefore, since e" 1 ^^ is bounded above and below on dB(xi, 1), we also get 

\K,rM\b(^i))\ <ce P 2 . 

We also require an estimate for the normal derivative of w l on dB(xi, 1). To get this, we use 
the slightly different representation 

J r J s 

Differentiating this gives 

\d r wi(x)\ < «!+ ^l))" 1 [ ^X(t)t N ^\Co(t)\dt for all x G dB(x u 1), 

and so 

sup \d rt w l (x)\ < cep 2 . (73) 

dB(xi,l) 

For j = 1, . . . , N, we may perform a similar analysis. Letting w l - denote the component of the 
solution in the j tn cigencomponent, then as before 

w){r) = Ki uRi {C mXi ,i))^Htt{ x - X i)/ R i) + 

where, denoting by 

= r ± ^ L (^^v ei (-log(r/R i )) - ^(-logW))) = log(r/i? 4 )), 

we represent 

J r J s 

The estimates 

\w)\<Cep\x~x t \ and |J% ifl< (C|B(* ( ,i))l < Cep (74) 

follow as above. 

The estimate for the normal derivative follows from the alternate representation 

J r J s 

which yields 

i^Mi < (^(i))- 1 [\l:%(t)t N -%(t)\dt. 

j Pi 

The estimate 

|5 n wj(x)| < cep^ 1 , Vx e aB(x 4 , 1), (75) 

follows immediately. 
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Let C 1 denote the sum over all higher eigencomponents of £. From (|70|) and (|7l|), and from 
Proposition |[ for each p, 6 (— iV, 2), 

l|G £i ,fl i (C i )l|2,a, M <Cep 2 -^ (76) 



where C is independent of e. In particular, using this with p = v, and also using ( |72| ) and (|74|), 

we get 

Finally, take p = in (76) to get the estimate 

\\d rj G^ Rl {C, l )\\c^(dB {Xt .i)) <cep 2 (77) 

on the boundary of B(xi, 1). 

Step 4. The final step is to consider the size of the correction term Wker, which is chosen 
to make the Dirichlet and Neumann data match. Its size is regulated by the size of the jump 
of the Neumann data of the interior and exterior solutions. Specifically, its Dirichlet data if? = 
(ipi, . . . , ip n ) 6 ®^ 1 C 1 ' a {dB{x i , 1)) is chosen to satisfy 

{ S e,R ~ %,r){^) = {(9 ri Wint,i - d ri w ext )\ ri=1 , (d rn w int , n - d Tn w ext ) \ rn=1 }■ 
Thus, we can estimate the size of if? from (|73|), (|75|) and (|7^), and the fact that the norm of the 



right side is bounded by Cep 2 . Hence ||ttffe e r ||2,o:,i/,2 < Cep as well. The estimate < Cep 2 

then follows easily. □ 

Proof of Lemma |^: The details of this proof are very similar to those of the proof above. 
In fact, proceeding through the steps of the previous proof, we find first that w ex t = 0. The 
solutions Wi n t^i may be found by the same procedure. The integral representations of Wj show 
that for j = and j = 1, . . . , N, the C 2 ,^ norms of these terms may be bounded by p 2 ~ u sup \ h Q \ 
and p 2 ~ v \ sup/i'l, respectively. The remaining details are left to the reader. □ 



10 The nonlinear fixed point argument 

We are now in a position to find a solution of the nonlinear problem. As usual, we fix the singular 
set E = {.Ti, . . . , x n } and positive parameters qi, . . . ,q n > associated to a set of admissible 
Delaunay parameters s. We also fix an approximate solution Ue{R 1 a, x) associated to the data R 
and a satisfying ( |l2| ) and (|l3|). The perturbation of u s (R,a,x) to an exact solution will involve 
not only a decaying term, but also a slight adjustment of the parameters R and a. Thus, recalling 
that an element w in M. has components (v, S,o), and changing our previous notation slightly, 
we wish to solve the equation Af(w) = where by definition 

N(N -2) - - n+2 
M(w) = A(u £ -(R + S,a + a 1 -)+v) + — K — ! -{u s {R + S, a + a, •) + v)~ . (78) 



Recalling the linearization A g ^ a of J\f at w = 0, we can rewrite (78), using a Taylor expansion, 
as ^ 

7VH = JV(0) + L Rra w + f (I>JVj( tOl s, a )-^(0 > S,a))(«.0,0) < ft 

Jo 

+(DAf l{0iSta) - DjV m o t o))(v, 0, 0) + f (DjV mStts) - DJV\ (0 , , ))(0, S, a) dt. 

Jo 

We have used a somewhat more elaborate expression than usual for the remainder term; this is 
necessary in the estimates below. Therefore, the equation J\f(w) = can be written as 



w 



A J,, AM): A (J (DjV KUl>s , s) - D/V^g^Xv, 0, 0) dt 
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+(#A/](o,s,a) - DAf li0 , 0fi) )(v, 0,0) +^ (D-Vj ( o,ts,t8) -DM m0t0) )(p,S,a)dtj 
We define the mapping V : M — ► M. by 

7» = -A-t _AT(0) - A-t _ (J (Dtf Kt0 ,S,*) ~DAf mStS) )(v,0,0)dt 



+ (DM m s,a) ~ DM K0>0fl) )(v, 0,0) + / (rWj (0 , t $, ta) -£>^j (0 .o,o))(0,5,a)dt 

Jo 

We shall prove that for e small enough, V is a contraction on a small ball in Ai of radius Coep 2 , 
where Co is twice the constant defined in Lemma [w]. This will give the existence of a solution of 
Af(w) = in the space M. 

Let u>i, w 2 be two elements in this ball. We begin by estimating the image by AT^ _ of 

I(w u w 2 ) = (lWj(« A , aa ) ~i^(o,3^ 

Lemma 11 There exists a constant c > (depending on Cq) such that 

ll^kr^'^IUf - C P 2 ~1K ~ w i\\m- 
Proof : From the structure of Af(w) we see that 

N(N + 2) 



((ug(R + S, a + a, ■) + v) "-^ - (u s (R + S, a + a, ■)) N ~ 2 )v. 
I(w 1 ,w 2 )(x) = N{N A +2 \ h+h+h), 



4 

We write 

N{N + 2) 



where 

h = ((ue(R+ S 2l a + a 2 ,x) +v 2 (x))t^ - (u e (R + S 2 , a + a 2 , x))^)(v 2 - Vi)(x) 
4 



N-2 



l 



v 2 (x)(v 2 — v±)(x) / (u g (R + S 2 , a + a 2 , x) + sv 2 (x)) N ~ 2 ds, 



(i 



I 2 = {{u g (R + S 2l a + a 2 ,x) + u 2 ) N - 2 - (u e (R + Si, a + a x ,x) + v t ) N - 2 )vi(x) 
vi(x) / [(ug(R+ sS 2 + (1 — s)Si, a + sa 2 + (1 — s)ai,x) + sv 2 (x) + (1 — s)vi(x)) ~" 



N-2 



N-2 







X 



22fj,i(R + sS 2 + (1 - s)S!,a + sa 2 + (1 - s)a>i, x)(S 2 ,i - Sx,i) 

i 

+ (R + sS 2 + (1 — s)Si,a + sa 2 + (1 - s)ai,x)(a 2i - a{ ti ) + s(v 2 (x) — vi(x))) ds 



*,3 

and finally 



I 3 = -((ug(R + S 2 ,a + a 2 ,x))"-i - (u g (R + Si, a + o>i,x)) v-*)vi(x) 
■-— ^— f ((u g (R + sS 2 + {l-s)Si,a + sa 2 + (l-s)cti,x))^ 



N-2 

(} jii(R+ sS 2 + (1 - s)Si,a + sa 2 + (1 - s)«i, x)(S 2 ,i ~ 
-^7i (-R + + (1 - s)Si,a + sa 2 + (1 - s)ai, x)(a^ - a^)) ds- 
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• First consider the region where dist(x, E) < e N - 2 . Here we use the fact that for some 
constant c > depending only on N, we have the estimates 

\fii(R+S, a+a)\ < cug(R+S,a+a) and \j?(R+S, a+a)\ < cdist(x, Ti)ug(R+S, a+a), 

where the constant c > depends on Cq. These yield 

\h\ <csp 2 -"snp(l,e^)\\v 2 - ui|| 2 , Q , 1/ , 2 dist(a;,S)^ +21 ' 
\h\ < cep 2 -"(dist{x,zy- 2 (\S 2 - Sx\ + dist(x,E)|a 2 -ai|) 
+ sup(l,e^)dist(a;,I]) i V !i+2,y ||u2 - Vi\\ 2 , a , v ,2) 

and 

141 < c£ / 9 2_l/ (dist(a;,I])' y " 2 (|52 - S x \ + dist(z, E)|a 2 -a x |). 

• Next, assume that dist(a;, E) e [e"- ! , 1]. In this case we use the estimates 

\m(R + S,a + a)\ < cedist(a;,E) 2 ~ A ' and \^ (R + S , a + a)\ 

< cedist(x, E)l + e^dist(a:, E) _Ar ), 
where the constant c > depends on Co- Then as above, we derive 

< ce^p 2 ~ v \\v 2 - ^H^^dist^E)^^, 

I Jii| < cep 2 - u (E^dist{x,i:) u - A (\S 2 - Si\ + (dist(a, E)"" 1 + e^dist^, E)" 1 )!^ - fii|) 

+e^dist(x, ^) N ~ 6+2v \ \v 2 - vi \\ 2 , a , v , 2 ) 

and 

|J 3 | < ce^p 2 - l/ (dist(a;,E) l/ - 4 (|52 - §i\ + (dist(x, E)^ 1 + £™=2 dist(a;, E) _1 )|a 2 — a\\). 

• Finally, if dist(x, E) > 1, we obtain 

|7i| <ce^C p 2 ~ v \\v 2 -i;i|| 2 ,a, v ,2dist(a;,E)- JV - 2 , 
\h\ < cep 2 -»(e^\S 2 - S x \ + e^\\v 2 - Vl \\ 2 , a ^ 2 )dist(x,i:)- N - 2 

and 

|/ 3 | < ce^C p 2 ^\S 2 - 5i|dist(a?, E)-^ 2 . 

Combining these gives 

||7(wi,w 2 )||o,a^-2,-2 < cp 2 ' u (e^\\v 2 - Ui | \ 2 , a<v , 2 + e(\S 2 - Si I +e^|a 2 -ai|)). 
In particular 

||.T(Wl, W 2 )\\ 0>a ,v- 2 ,- 2 < Cp 2 ~"(jf\\v2 ~ Ul||2,a,i/,2 + e|^2 - <§l| + £p|a 2 ~ <5l|) 

= cp 2 ^||w 2 -u;i||_^. 

The estimate of the Lemma now follows from the boundedness of AT^_, as proved in Proposi- 
tion!. ' n 

Next shall estimate the image by A~]^ of 

//K,io 2 ) = (DM m§2 ^ 2) -DM mo , Q) )(v2,0,0) - (DJVj (0 ,s 1)Sl ) -^(o,o,o))(^, 0,0). (79) 
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Lemma 12 There exists some constant c > (depending on C ) such that 

I \Kk« H ( Wl ' ^) I U < CP 2 '" | \w 2 - Wl | \ M . 

Proof : Recall that 

^(^ + 2) - - 4 4 

{DM mSA) - DM mm ){v, 0, 0) - + 5, a + a, •))*=* - (u £ (fl, a, ■))*=* )«• 



Therefore, we may decompose II as 



where 



Ih = ((ue(R + S 2 ,a + a 2 ,x))"-2 - (u s (R,a,x))tr=s)(v2 - vi)(x) 



and 



4 f 1 _ - 6~iv 

= jj— ^{v 2 - vi){x) J {{u g {R+ sS 2 ,a + sa 2 ,x))»-* 

(Y^ Hi{R + sS 2 ,a + sa 2 , x)S 2ji + ^ ll{R + s $2, a + sa 2l x)a{ A ) ds) 

i i,j 

Ih = {(ue(R+ S 2 , a + a 2l x))^2 _ (u s {R + Si, a + an, x))~^)vi(x) 

4 f 1 _ - - 6-A 

= N _ 2 vi(x) J ((iie(R + sS 2 + (l-s)Si,a + sa 2 + (l-s)ai,x))~^ 

, a + sa 2 + (1 — s)ai,x){S 2 A ~~ 

+ Xl 7 i ^ + S ^ 2 + s )Si,a + sa 2 + {1- s)a 1 ,x)(ai ti -a{ ti ))ds. 

• As before, first consider the region where dist(x, S) < e N - 2 . As in the previous Lemma, 
we use the fact that there exists a constant c > only depending on N such that 

\p>i(R + S, a + a)\ < cu g (R + S, a + a) and < cdist(x, T,)u s (R + S, a + a). 

So, we get the estimate 

|/7i| < cdistfo T,) u - 2 (C oP 2 + dist(x, T,)C p)\\v 2 - «i||o,a,.,,2, 

\II 2 \ < cep 2 -" d\si{x, Y,y- 2 {\S 2 -Si\ + dist(:r, £)|a 2 -ai|). 

• Next, assume that dist(x, S) e [s 7 ^ , 1], then as above, we can derive the estimate 

|//i| < ce^dist(a;,E) ,y - 4 ( / 9 2 + (dist(a;, E)^ 1 + e^dist(a;, S)" 1 ) / 9)||w 2 - «i||o,a, v ,2, 
|// 2 | < cp 2 - ,y e^dist(a;,S) ,y - 4 (|5 2 + (dist(x, E)^ 1 + e™^2dist(x, E) _1 )|a 2 — ot\\). 

• And finally, if dist(x, £) > 1, we obtain 

|7Ji| < ce^dist(a;,S)- Ar - 2 / 9 2 ||w 2 - wi|| , Qjl ,, 2 , 
\U 2 \ < cC Q p 2 - v e^ dist(x, Y,)~ N ~ 2 \S 2 - 5i|. 
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This allows us to conclude that 

\\II{wi,w 2 )\\o, a ,v-2-2 < cC Q (p 2 \\v 2 -v 1 \\2, a ,v,2+p 2 ~ v e(\S2-Si\ + e^\a 2 -ai|)), 
which leads to 

\\II(wi,W 2 )\\o,a,v-2,-2 < cC p 2 ~ U \\w 2 ~ WlWjft. 

Then the estimate of the Lemma follows from the result of Proposition |j| □ 
We now turn to the estimate of the image by At^ _ of 

III{ Wl ,w 2 ) = (rW[( 0l s aj a a ) -^(o,o,o)) (0,5 2 ,a a ) - (Dtf m § 1>Sil) - DJVj (0) o,o))(0,Si,ai) 
As before we decompose III = III\ + III2 where 

Ilh = {DM ms ,. &2) - DN mm ){Q,S 2 -Suaa-at) 

Ilh = (DM mS ^ 2) - DX mSl , &1 ))(0,Si,ai) 

To estimate either of these terms, it is convenient to consider them as special cases of the more 
general quantity 

III(S, S',S",a,a',a") = (DJV^ 0< §' +5, - D N\(o, S", a") 
We first prove the Lemma : 

Lemma 13 There exists some constant c > ( depending on Cq ) such that 
Proof : We may write 

»i 

5' an ~> ^.ii\ I 7-1 2 , 



III(S, S', S", a, a', a") = / D Mr 0j s+ s s>,<x+sa')(( ' 5 "' a '), (0; S", a")) ds, 

Jo 

and shall concentrate on estimating the integrand 

iit(S, S', S", a, a', a") = £> 2 A/j (o ^ a) ((0, S', a'), (0, S", a")). (80) 
We take advantage from the observation that, for all R 6 K™ and for all a G (SSL N ) n 

N(N-2) - n+2 
Au & (R, a, x) H x — J -(u g (R,a,x))^ = 0, 

if dist(x, Xi) < pi. Therefore, if dist(z, ac<) < pi, (^) is identically 0. Moreover, if dist (x, Xi) > 2pi, 
we observe that Ue(R, a, x) = Wg[R, x) is harmonic and does not depend on a. Therefore, in this 
case, Hi reduces to 

m = Y, d k, Rs + s) «s )sis>;, 

And, still assuming that dist(a;, Xi) > 2pi, we find the estimate 

\iii\ < ce^d\st(x,-E)- N - 2 \S'\\S"\. 
Finally, if dist(a;,ari) € [pi,2pi], then by definition 

u g (R,a,x) = Yxi( x - Xi)u Si (Ri,ai,x- x t ) + wg(R,x)(l - Y^Xijx - Xi)), 
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so we can write 

. , - _ s N(N-2),_ ,- _ .. iv+2 
Au g (R,a,x) + — ^ ^(u £ -(i?,a,x))«-2 

N(N — 2) W + 2 _ - iV + 2 

= — ^ (Ri,a,i, ■ ~ Xi ))^ + (XiU £i (Ri, o,-, • - x 4 ) + (1 - x»)fr(-R, O) 1 ^) 

+A(xi)(u £i (Ei,ai, • -Xi) -We(R,-)) + 2V(xi)V(u Ei (-Ri, a», • - a;,) -u)e(E, •)). 
Differentiating twice with respect to both _R and a, it is easy to derive the estimate 

l^| ( S+s 1 a+a)(-Xi(« s «)^ + (XiU Si + (l-Xi)w s )^)(S',a')(S",a")\ 

< ce^(p-\p 2 - N \S'\\S"\ + p(\a'\\S"\ + \a"\\S'\) + p 2 \a'\\a"\) 
+p N - 6 (p 2 - N \S>\+ p\a>\)(p 2 ~ N \S"\+ p\a"\)) 

< ce^(p- 2 - N \S'\\S"\ + p- 3 (|a'||5"| + \a"\\S'\) + p' 2 \a'\\a"\). 

Using the result of Proposition [| and differentiating twice with respect to both R and a, we get 
the expansion 

T/ \(R+S,. 

= D 2 iR+Sra+ ^( S fR^f £ f -x t \ 2 - N )(S>,a>)(S",a") 

AT _ O 2 — TV p . JV-2 _ _ 

^ + 5, 5+a) (feo^ 2 s + Ei fl . 2 k^.r N k^ i ))-(^-^))(5' I «')(s''a") 

+ 0( £ p 2 )(|5'||5"| + |a'||a"| + \S'\\a"\ + \a'\\S"\). (81) 

and 

D \{R+S,&+&)^ dr *a ( U e> i R ioi a io> x ~ x io) ~ Wg(R, x)))(S' , 5c')(S" , 6l") 



D " B < s,a+a)( u ^ -w g (R,x))(S',a')(S",a") 



N — 2 2-N ^ £^ N-2 X — Xi — f - f/ , f 

= ~^r- D \(R+5,a+ci)(( e io R i 2 a io-2_^T R * 2 \ X io~ x i\~ (^o ~ x i)) ' [~ _ *° i ) (S', « « 

+ 0( e p)(|S'||S"| + \a'\\a"\ + \S'\\a"\ + \a'\\S"\). (82) 
This expansion allows us to estimate 

'\(B.+ S.a+s)( A (Xi)(Ue,(Rt, CLi, ■ - Xi) - W S (R, ■)) + 

2V( Xl )V(u £t (Ri,ai, --Xi)- w g (R, ■)))(§', a')(S", a"). 



D 2 



As in the setting of Lemma ^, we shall decompose this quantity as h + h' + h" as before, then 
we get the estimates 

\ho\ < ce{p~ 2 \S'\\S"\ + \a'\\a"\ + \S'\\a"\ + \a'\\S"\). 
\h'\ < ce(p- 1 \S'\\S"\ + \a'\\a"\+ p- 1 \S'\\a"\+ p-^a'llS"]). 

and finally 

\h"\ < ce(\S'\\S"\ + \a'\\a"\ + \S'\\a"\ + \a'\\S"\). 
Therefore, using the result of Proposition ^ as well as the result of Lemma ^| we obtain 

\\A~ l &s iii(S,S',S'',a,a',a")\\j^ < ce-(|5'||5"| + p 2 \a\\a"\ + p\S'\\a"\ + p\a'\\S"\). 
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As a consequence we get that 

□ 

As a corollary of this last lemma, recalling that e|5j| + £p|<5j < C^ep 2 , we get that 
WKka IH ( Wl ' W2 ^M ^ C£ (P 2 |^2- Si| + p 3 |a 2 -ai|) < cp 2 ||w 2 - ioxII^. 

Collecting the results of these lemmata, we have established that for some constant c > 
depending on Co, 

\\V(w 2 ) - < cp 2 - v \\w 2 - Wl \\ M . 

By taking e, and hence /?, sufhciently small, we have shown that the map V is a contraction on the 
ball of radius C^ep 2 in M.. Hence it has a unique fixed point w, and this function w — (v, S, a) is 
a solution of the equation (|78|). It is clear that Ug(R 4- S, a + a, •) + v is positive near the points 
of S, hence by the maximum principle is positive everywhere. 

This completes the existence of the solution promised in Theorem [p. 

11 The nondegeneracy of the solutions 

We now show that for e sufficiently small, the solutions we have constructed above are non- 
degenerate in the sense defined in 0, Q. Actually, there are two closely related notions of 
nondegeneracy, the definitions of which we now recall, in terms of the notations of this paper: 

4 

Definition 3 Let go be the standard metric on S . A metric g = u N ~ 2 go of constant positive 
scalar curvature on S \ A, as well as the corresponding conformal factor u, is called marked 
nondegenerate if the linearization of the scalar curvature operator Af for go about the solution 
u is injective on the function space C% a (S N \ A) for all fi > (2 - N)/2, or equivalently, if the 
linearization of M relative to the metric g about the constant solution 1 is injective on the function 
space C^, a (S N \ A) (defined with respect to the metric g) for all // > 0. The metric g, or solution 
u, is called unmarked nondegenerate if the linearization is injective for all ji > (4 — N)/2, or 
equivalently, for all fi' > 1. 

These two nondegeneracy conditions are precisely what is needed to ensure the smoothness 
of the marked and unmarked moduli spaces A4\ and A4 n at g. The former of these spaces is 
the space of all metrics of constant positive scalar curvature on the complement of the finite 
set A in the sphere S N , while the latter is the set of all such metrics on the complement of 
any finite set A of cardinality n. As proved in [Q and [§j, these moduli spaces are real analytic 
sets, hence are stratified and may be written as the union of smooth, real analytic manifolds 
of varying dimensions. The existence of one smooth point in a given component shows that 
the top dimensional stratum in that component is of the dimension predicted by the formal, 
index-theoretic, calculations, namely n for the marked spaces and n(N + 1) for the unmarked 
ones. 

In this section we shall prove 

Proposition 10 The solutions constructed here are unmarked nondegenerate. For generic con- 
figurations A they are also marked nondegenerate. 

As a corollary, we obtain smoothness of the unmarked moduli spaces without any restriction 
(when e is small) and of the marked moduli spaces for generic configurations. In fact, near 
a generic A, and when e is sufficiently small, we may use (pi, . . . ,p n , ei, . . . , e n ) — {p, e) € 
(S N ) n x (K+) n as coordinates on the unmarked moduli space M n , and e as coordinates on the 
marked moduli spaces. We note that marked nondegenerate solutions (which are a fortiori also 
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unmarked nondegenerate) were constructed in || for certain very special configurations A, which 
in particular contain only even numbers of points; the Delaunay parameters of those solutions 
are not precisely prescribed, but they need not be close to zero. Even when n is even, it is not 
clear that these solutions lie in the same component of M n as the ones we construct above. 

We first demonstrate the unmarked nondegeneracy; the proof is by contradiction. By a slight 
change of notation from the rest of the paper, we consider the solution u on S N \ A to have the 
form u — a, •) + v, where the approximate solution on the sphere (rather than on M. N ) is 

now denoted v,g(R,d, •), and v is an element of C 2,Q with 1 < v < 2, with norm in this space 
bounded by C§ep 1 ~ lJ . Assume that for some sequence of tending to 0, the linearized operator 

N(N — 2) , N(N + 2) ^ 
L k = A s n + u k 

is not injective on C^' a (S N \A) for some [i > (4 — N)/2. Here we have denoted by Uk the solution 
Ug k (Rk 7 a/c, •) + Wfc. Then there is some element Wk € C^' a such that CkWk = 0. 

First normalize Wk, multiplying it by a suitable constant, so that sup d(y)~ tJ '\wk{y)\ = 1, where 
d(y) is the distance of the point y from A, say in the spherical metric. Choose a point yk 6 S N \A 
realizing this suprcmum, i.e. such that c?(j/fc ) — M |iffc (yfc ) | = 1- As k tends to infinity, the function 
Vk tends to zero in C^' Q , since its norm is dominated by Co£kpl~ u , an d % converges uniformly 
to zero on compact subsets of S N \ A. In addition, we can assume that the sequences Rk and dk 
converge to some fixed R* and a* , respectively. 

Suppose first that some subsequence of the yk converges to a point yo G S N \ A. Then we 
can extract a subsequence of the Wk which converge (in C°°) to a limiting function w on S N \ A; 
w must be nontrivial since d(yo)~ tJ, \w(yo) = 1. Furthermore, Lw = (A s n — (N — 2) 2 /4) w = 0, 
and \w(y)\ < d(y)^ . Since fx > (4 — N)/2, it is standard that w is a weak solution of Lw = 
on all of S N , hence also a strong solution. But clearly L has only trivial nullspace, hence w = 0, 
which is a contradiction. 

If, on the other hand, some subsequence of the yk converge to one of the points Pi £ A, then 
we must argue somewhat differently. Choose a function A on S N \ {pi , <?}, where q ^ A, which 
transforms this twice-punctured sphere to the cylinder K x S 1 ^^ 1 and such that the standard 
spherical metric gs and the (dilated) cylindrical metric gc are related by 

g s = A7^g c , g c = ~^{dt 2 + d9 2 ). 

2 JV Tir 2 — N 

On C, the function A is simply a multiple of (coshi) - , and on S is of the order d(y,p)^~ , 
for p — pi or q. Noting that gc has scalar curvature N(N — 1), the same as gs, the conformal 
Laplacians of these two metrics, Lc and Ls satisfy the usual transformation rule 

L c (Adp) = A^L S (j>, (83) 

for any function <f>. The solution Uk on S corresponds to a solution Auk on C. It is easy to 
check that the linearizations £s,k of the scalar curvature operators on S at uu and £c,k on C at 
Auk satisfy the same transformation rule ( |83] ) as is and Lc-, cf. (Unlike the transformation 
rule for the conformal Laplacian, this holds only because we are dealing with two metrics of the 
same constant scalar curvature.) 

Because of these machinations, we may replace Wk by a solution Wk of Cc,kWk = 0. For 
convenience, we call Cc,k simply Ck and Wk simply Wk again, and let y — (t, 9) denote the 
variable on the cylinder. Define p! = fi + (N — 2)/2, so that // > 1. These new functions Wk 
satisfy 

sup%)^V fc | = l, (84) 

where d{y) is some function equalling the distance to the set of other singular points A' = A\{pi } 
(transplanted to the cylinder) in a neighbourhood of this set, and equalling secht outside this 
neighbourhood. Because of (J8J), Wk decays at either end of the cylinder. 
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As before, let — (tk,&k) denote the point on C where the supremum in (84) is attained. 
We already are assuming that tk — > oo. By translating back by tf. and multiplying by a suitable 
constant, we find yet another sequence of solutions, which we again call Wk, attaining their 
maximum at t = 0, and which solve the translated equation, which we again write as CkWk = 0. 
Here Ck is the linearized scalar curvature operator relative to a metric g^ — (A(ug k +Vk)) 4 ^ ( N ~ 2 ^ gc 
which is singular at a finite collection of points which are translating toward t = — oo. Since Vk is 
tending to zero, gk more and more nearly approximates the Delaunay metric Ve k ^ 9c, where 
V €k is one of the (0-independent) Delaunay solutions considered in §2. 

It is not hard to see that the Wk converge to a nontrivial solution w of the limiting equation 
Cw = 0, and that w is bounded by e _M * for all t. 

There are two subcases. In the first, T4(0, 6*) tends to zero. Then w satisfies the equation 

N A \ N(N-2) N / , (N-2) 2 \ 

JT2 ( 9 * + A ^"0 w ~ 4 w = W^2 \ di + ~i) W = °- 

By decomposing w into eigencomponents with respect to A s n-i , we see that any eigencomponent 
Wj is a sum of exponentials, Wj — a+e 7jt + aje~ ljt . Since w decays as t — > +oo, a+ = 0. But 

then it is clear that no function of the form e _7j t can be bounded for all t by e _/i * unless \j! = jj, 
which is not the case, so we arrive at a contradiction. In the second, 14(0, 9) does not tend to 
zero. By translating by a fixed finite amount, we may assume that Vk tends to the function 
(coshi)^ 2- ^/ 2 , and hence, after pulling out the superfluous constants, that the limiting function 
w satisfies 

n9 * (N -2) 2 N 2 -A 2 \ 

d 2 + A S n-i - - — + — - — scch 2 t I w = 0. 



n2 f(N-2) 2 \ N 2 -4 2 

d t w i - a + X i W 'J + — a — sech 1 w j = °- 



Again separate w into its eigencomponents wj . Then 

For j = the indicial roots of this equation at both ±oo are ±(N — 2)/2, for j = 1, . . . ,N they 
are ±N/2, and for j > N they all satisfy (t^I > (N + 2)/2. 

The components Wj with j > N are easy to eliminate. In fact, these Wj must decay faster 
than e ± ^ JV+2 -" t '/ 2 at ±oo, so we may multiply the equation satisfied by Wj and integrate by parts 
to obtain 

/ (N — 2) 2 \ N 2 — 4 
{dtWjf + iXj + \ w j 4 — sech 2 i w 2 dt = 0. 

Since Xj > 27V, the integrand is nonnegative, hence Wj = 0. 

For the remaining cases, when j < N, the indicial roots at ±oo are less than (N + 2)/2 in 
absolute value. On the other hand, to check unmarked nondegeneracy it suffices to use any fi' > 1 
which is also less than the next omitted weight for the linearization about a Delaunay solution, as 
determined in §5. But by the results of that section, as e tends to zero, this next omitted weight 
tends to (N + 2)/2. Thus now choose p! in the range (N/2, (N + 2)/2); then each of the Wj with 
j < N decays less quickly than * as t — > +oo, which implies that these Wj too must vanish. 
This is a contradiction. 

This covers all cases, so we have showed that the linearization is injective on the appropriate 
weighted Holder spaces. 

This completes the proof of unmarked nondegeneracy of the solutions when e is sufficiently 
small. We complete the proof of the proposition by demonstrating the marked nondegeneracy of 
solutions for generic configurations A. This is a simple consequence of Sard's theorem. In fact, as 
discussed in j^], near smooth points of the unmarked moduli space Ai n , there is a real analytic 
fibration 7r : M n — > C n onto the configuration space of n distinct points in S . C n is itself a 
real analytic manifold. Standard differential topological arguments now imply the surjectivity of 
the differential 7T* for all points in a generic fibre of tt. By a dimension count, surjectivity of this 
differential is equivalent to the marked nondegeneracy of all points in the fibre. 
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12 The general case 



In this brief final section we discuss the essentially minor changes that need to be made in order 
to prove the more general statement of Theorem [l], where the singular set is allowed to have 
components of positive dimension. As in the main body of the paper, there are three steps. 
First, we must construct an approximate solution, or rather, a family of approximate solutions 
that become increasingly concentrated at A; next we must prove that the linearization about 
one of these approximate solutions is surjective on an appropriate function space provided the 
approximate solutions are sufficiently concentrated; finally, we perturb once again to an exact 
solution. In each of these steps, we must somehow combine the constructions and proofs from 
our previous paper || with the ones here. The last step, showing that an appropriate map is a 
contraction, is straightforward and we shall not comment on it further. We now describe the first 
two steps. 

As in the introduction, divide the components of the singular set A into two groups, A = 
A'UA", where the first is a finite set and the latter contains all the higher dimensional components. 
The construction of a family of approximate solutions around a component Aj of dimension k, 
where < k < (N — 2)/2, is given in detail in ||, but briefly, it is obtained by first fixing a 
tubular neighbourhood T(Aj) and identifying it with a neighbourhood of the zero section in the 
normal bundle NAj . On each of the fibres N p Aj we glue in a sufficiently dilated solution of the 
equation 

N+2 

Au + u N - 2 = 0, 

cut off to be supported in T(Aj). Because the fibres are of dimension N — k, this equation is 
subcritical, and one can show that the radial solutions have the form \x\ < ' 2 ~ N ^ 2 v, where v tends 
to zero rather quickly as \x\ — > oo. There is a one-parameter family u £ of such solutions, which 
are all, in a suitable sense, dilated versions of one another. Again, we refer to |(| for all details. 

To construct the approximate solutions when A' and A" are both nonempty, we fix Delaunay 
parameters e' for the points of A' and dilation parameters e" for the components of A", and set 
e — (e',e"). The elements of each of these subsets are mutually commensurable, i.e. s' = e'q' 
and e" = e"q", where q = (q',q") is a vector with all components positive. Now construct the 
approximate solutions singular at the points of A' as in §3, balanced exactly as before (here we 
use that |A'| > 2). This approximate solution is of size e' outside a fixed neighbourhood of A', 
and so we use a partition of unity to glue it to the approximate solution defined in each TAj . 

The main step is to show that the linearization of the scalar curvature operator, L Sl is surjec- 
tive as a map 

U ■■ C%> (S N \ A) © W ► Cy_ 2y/ _ 2 (S N \ A), 

at least when all the components of e are sufficiently small. Here 1 < v' < 2 is the weight 
parameter determining growth of functions in a neighbourhood of each point of A' and v" < 
(2 — N)/2 determines the growth near each component of A". W is the same deficiency space as 
before. From the results of this paper, there exists an inverse for the linearization about the part 
of the approximate solution which is singular only at A'; we denote this inverse, as well as its 
Schwartz kernel, by H' g , . Let H",, denote the right inverse, or Schwartz kernel, for the linearization 
about the part of the approximate solution which is singular only at A", as constructed in |^. 
Now let H g be a Schwartz kernel obtained by using a partition of unity (e.g. the same one as 
used to construct the full approximate solution) to glue together these two pieces. Then it is easy 
to check that 

L g H g = I + R, 

where R is of size max{e',e"}. Clearly, then, / + R is invertible, and a right inverse for L g is 
given by H E (I + R) . It is straightforward to verify the necessary mapping properties for this 
operator. 

As noted above, the final step, using a contraction mapping argument to show the existence 
of a solution of the problem, is standard. 
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This completes the proof of Theorem |lj in full generality. 
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